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Apresentação 


Esta obra tem por finalidade trazer ao aluno o primeiro contato com sistemas de controle, 
apresentando conhecimentos elementares de matemática aplicada e de álgebra. O livro busca servir 
de apoio a outros materiais, procurando desmistificar a tarefa de se projetar um sistema de controle 
que opere sobre um dado processo ou planta. São fornecidos exemplos de aplicação de sistemas de 
controle, não se perdendo de vista a necessidade de o aluno buscar aprofundar sempre seus conheci- 
mentos matemáticos, que servem de forte embasamento à compreensão dos conceitos aqui apresen- 
tados, ao mesmo tempo em que se procura trazer um foco de compreensão fundamentado em regras 
específicas de projeto, claras e imediatas. 


Esta obra está dividida em onze capítulos. O Capítulo 1 versa sobre o histórico da Teoria de 
Controle, apresentando conceitos fundamentais desta área de conhecimento. O segundo capítu- 
lo apresenta o referencial teórico de manipulação de vetores e matrizes, úteis no entendimento de 
como se processam algumas das técnicas de controle. O Capítulo 3 traz as representações matemáti- 
cas mais usuais de sistemas, apresentando o ferramental necessário à sua manipulação, nos domíni- 
os do tempo e da frequência. O Capítulo 4 discute o aspecto da estabilidade de sistemas, focan- 
do-se naqueles de primeira e segunda ordens, posto ser esse aspecto a meta fundamental em todo 
e qualquer projeto de controle. O Capítulo 5 discorre sobre a técnica de Lugar das Raízes, relevante 
na caracterização da resposta de sistemas. Os Capítulos 6 e 7, apresentam, respectivamente, técnicas 
rápidas de projeto de sistemas de controle, baseadas na teoria clássica e na construção de Lugar das 
Raízes. O Capítulo 8 discute o projeto de controladores no domínio da frequência, com ênfase no 
uso de Diagramas de Bode. O Capítulo 9 apresenta a sintonia de controladores PID por meio das 
regras de Ziegler-Nichols. O Capítulo 10 apresenta os primeiros elementos de análise de sistemas de 
controle no espaço de estados, com uma introdução preliminar a projetos. O Capítulo 11, por fim, 
traz considerações sobre o uso de atuadores e sensores em sistemas de controle aplicados a plantas 
industriais. 


Os apêndices da obra trazem elementos matemáticos importantes ao entendimento da teoria 
apresentada. 


O autor 
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Introdução a 
Sistemas de 
Controle 


Para começar 


Este capítulo tem por objetivo apresentar os conceitos básicos pertinentes a sistemas de controle. 
O capítulo busca abordar, de forma didática, o conceito de sistema, trazendo à luz as definições de sinal 
de entrada e resposta, visando a um futuro entendimento de como um sistema de controle opera sobre 
tais sistemas da natureza. Ao mesmo tempo, o conhecimento histórico exposto consegue trazer ao aluno 
a justificativa de uso das técnicas de controle, enfatizando a importância da Teoria de Controle em enge- 


nharia e tecnologia. 


As informações apresentadas neste capítulo constituem uma importante base para a compreensão 
gradual de todos os conceitos inerentes a sistemas de controle, a partir das definições abrangidas e da 
metodologia de apresentação utilizada nesta área de conhecimento. 


1.1 Histórico da teoria de controle 


A noção de “sistema de controle” remonta ao conceito de necessidade de se observar a resposta 
de uma determinada estrutura a um estímulo de entrada, a fim de que tal resposta seja conduzida a 
um valor considerado ideal. Tais preceitos se iniciaram na Idade Contemporânea com o estudo da 
homeostase por Norbert Wiener, especificamente em meados do século passado, advindo, todavia, 
de tempos mais remotos, quando o homem já se preocupava com a medição precisa do tempo, pas- 
sando pela Revolução Industrial, pelas grandes guerras e ganhando corpo com o início da chamada 
corrida espacial, em 1957. 


A Teoria de Controle tomou corpo, de fato, a partir da concepção de “controle automático” O 
primeiro controlador automático industrial de que se tem notícia se refere ao regulador centrífugo 
concebido em 1775 por James Watts, que intentava prover o controle de velocidade de máquinas a 
vapor. Em 1868, James Clerk Maxwell deu um passo relevante nas aplicações da Teoria de Controle, 
ao analisar matematicamente um sistema máquina-regulador realimentado. Durante a Primeira 
Guerra Mundial, Minorsky desenvolveu o servomecanismo, também baseado em realimentação, que 
visava à manutenção automática da rota de navios. Em 1948, o trabalho supracitado de Wiener, des- 
crevendo fenômenos neurológicos e os sistemas de controle preexistentes no corpo humano, trouxe 
considerável progresso no entendimento de sistemas mais complexos. Em 1932, H. Nyquist, funcio- 
nário da Bell Labs., elaborou a primeira teoria geral de controle automático com a chamada Teoria da 
Regeneração, pela qual se estabelecia um critério específico para o estudo da estabilidade de sistemas. 


É fato que muito do avanço da Teoria de Controle adveio justamente da compreensão cada vez 
mais apurada que o homem tinha de seu entorno, envolvendo as noções de medição de tempo e de 
espaço, bem como de variáveis de processo como temperatura, vazão, pressão, nível, entre outras. 
Do mesmo modo, a História apresenta relatos de desenvolvimento de reguladores e comparadores, 
sempre tomando como base o conceito de realimentação, conceito este que atingiu tecnologicamente 
a produção de máquinas, veículos e dispositivos criados pelo homem. Toda a discussão em torno 
do controle de sistemas tinha fonte, ao longo dos tempos, da própria definição de “sistema”, situada 
no domínio do conhecimento humano e referente a uma entidade dinâmica sujeita a excitações de 
entrada e respostas de saída, representativa de um fenômeno natural. 


Dentro do campo de análise de sistemas, diferentes estratégias também se fizeram surgir. Téc- 
nicas baseadas no domínio do tempo foram confrontadas com outras baseadas no domínio da fre- 
quência, que ganharam força com técnicas de controle destinadas à remoção de ruídos. Técnicas 
baseadas na definição de funções de transferência, diagramas de blocos e lugar geométrico das raí- 
zes ensejaram, assim, um avanço na caracterização matemática de respostas de sistemas a entradas 
conhecidas. Na passagem das técnicas de controle clássico para as de controle moderno, abandona- 
ram-se as técnicas no domínio da frequência e surgiram novamente algoritmos baseados no domí- 
nio do tempo. Igualmente, o advento de computadores digitais trouxe benefícios variados à Teoria 
de Controle Moderno, no sentido de que a resolução de equações associadas a problemas de controle 
se podia fazer de modo off-line. Além disso, características de variância no tempo podem ser mais 
bem tratadas mediante o emprego de técnicas de processamento digital de sinais. 


Fique de olho! 


O lançamento do satélite Sputnik pelos soviétivos, disparou uma verdadeira corrida norte-americana por soluções diferen- 
tes em controle de trajetória de veículos aeroespaciais, já que se quebraram paradigmas da Teoria de Controle Clássico, 
motivando um retorno às técnicas antigas de controle — vindo daí o retorno à modelagem de sistemas no domínio do 
tempo! 


Enfim, nos dias atuais, há uma conjunção de técnicas de controle clássico e moderno, de modo 
a sanar fragilidades deste último grupo ao não fazer uso da intuição proporcionada pelas técnicas 
baseadas no domínio da frequência, além de não apresentar a mesma robustez daquelas técnicas. 
Assim, fundem-se vantagens de ambas as escolas, convergindo para uma nova teoria de controle que 
agrega as melhores características de técnicas clássicas e modernas. 
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1.2 Noções teóricas básicas 


Sinal de Sinal de 
entrada saída 


Sistema 


Figura 1.1 - Representação gráfica de um sistema. 


1.2.1 Definição de sistema 


Define-se sistema como a representação de qualquer fenômeno natural, físico ou químico, 
sujeito a estímulos externos (ditos sinais de entrada ou variáveis de controle) e produzindo respostas 
a tais estímulos (ditos sinais de saída ou variáveis controladas), como pode ser visto na Figura 1.1. 
Um sistema, de modo peculiar, simboliza um conjunto de elementos intrinsecamente conectados, 
formando um conjunto internamente organizado. De modo geral, um sistema é dotado da capaci- 
dade de gerar uma saída específica a uma entrada conhecida (saída esta exata ou ao menos aproxi- 
mada), de maneira que a integração de suas partes conduz a um resultado combinado dos efeitos 
individuais de cada porção. São exemplos de sistema um foguete em sua trajetória espacial, um auto- 
móvel em movimento ou mesmo um chuveiro aberto, derramando água durante seu banho! 


Matematicamente, um sistema pode ser classificado como linear ou não linear. Diz-se linear 
aquele sistema regido por finitas relações matemáticas entre saída e entrada dispostas em também 
finitas equações lineares - em outras palavras, por equações algébricas que contêm apenas variáveis 
de primeiro grau, ou por equações diferenciais lineares obtidas de leis físicas que governam um dado 
sistema em particular, de modo que vale o chamado princípio da superposição (a resposta para 
vários sinais de entrada equivale à soma das respostas individuais a cada um dos sinais de entrada). 
Por sua vez, sistema não linear é aquele cujas variáveis não obedecem às condições de linearidade. 
Em outras palavras, um sistema é chamado de linear quando é representado por equações lineares: 
matematicamente, é aquele em que, se x, e x, são soluções de tal sistema e cj e c, são constantes 
arbitrárias, logo c;x, + c,x, também é solução do sistema - aplica-se o chamado princípio da super- 
posição. Por exemplo: 


2x, +3x, + 4x; =12=> compõe um sistema linear 


4x2 +3x,? + 5x3? =16 = não compõe um sistema linear 


Deste modo, a condição de linearidade de sistemas descreve uma categoria particular de siste- 
mas a que as técnicas de controle descritas nesta obra se aplicam. É possível, por meio da linearização 
de sistemas não lineares, objeto de seção descrita no Apêndice desta obra, proceder às mesmas técni- 
cas utilizadas para sistemas lineares. 


1.2.2 Termos usuais em sistemas de controle 


Além da definição formal de sistemas, apresentada na seção anterior, há um glossário de ter- 
mos que devem ser rigorosamente assimilados para que a compreensão acerca de sistemas de con- 
trole se faça mais evidente. São eles: 
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» 
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» 


» 


» 


» 


» 
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Planta: dispositivo ou equipamento que tem por finalidade executar dada tarefa (por 
exemplo, motor, braço robótico, entre outros). 


Processo: fenômeno, natural ou artificial, que obedece a leis matemáticas bem descritas. É, 
na verdade, a essência do funcionamento de um sistema (por exemplo, fermentação, des- 
tilação fracionada, queda de preço de uma mercadoria, entre outros). 


Especificações de desempenho: descrições do comportamento a ser apresentado pelo sis- 
tema, de acordo com o desejado pelo usuário. 


Sistema: conjunto ou coleção de partes, dentro de um universo maior, que estão interligadas 
ou relacionadas de tal modo a constituírem um todo. 


Controlador: dispositivo empregado para se perfazer o controle de um sistema, ou seja, 
atingir as especificações de desempenho desejadas. 


Sistema de controle: conjunto formado por sistema a ser controlado e controlador. 
Setpoint: valor desejado para o sinal de saída de um sistema, dada uma entrada conhecida. 


Malha aberta: categoria de sistema (ou de malha de controle) em que o sinal de saída 
depende apenas das características do sistema e dos sinais de entrada deste sistema (por 
exemplo, veículo sem pedal do freio, máquina de lavar). Em outras palavras, a ação de 
controle independe da saída apresentada pelo processo que se deseja controlar. 


Malha fechada: categoria de sistema (ou de malha de controle) em que o sinal de saída 
depende não somente das características do sistema e dos sinais de entrada deste sistema, 
como também dos sinais prévios de saída (por exemplo, veículos de carga com velocidade 
controlada por dispositivo). 


Distúrbio ou perturbação: todo e qualquer sinal que venha a provocar interferência ou dis- 
torção nos sinais presentes em um sistema (por exemplo, ruído acústico, “manchas” em 
imagens transmitidas, perdas de bits em pacotes de dados de uma rede, radiação eletromag- 
nética interferente em um sinal de rádio). Em uma malha de controle, são especialmente 
importantes as perturbações de processo (load disturbances) e as perturbações de setpoint. 


Controle preditivo (feedforward): estratégia de controle realizada unicamente com base nos 
dados de entrada. Para sua aplicação, o controlador deve entender as relações de causa e 
efeito relativas ao comportamento do processo. 


Realimentação (feedback): laço de “retorno” do sinal de saída à entrada, com o fim de ele- 
var (realimentação positiva) ou reduzir (realimentação negativa) o nível de sinal a alimen- 
tar a planta - e, portanto, buscar a obtenção de um sinal de saída desejado. Nessa estratégia 
de controle, é realizada uma comparação entre o resultado obtido e o desejado. 


Ordem de um sistema: ordem da equação diferencial que o representa. 


Os conceitos de malha aberta e malha fechada devem ser detalhadamente descritos e explora- 
dos, sendo objetos das seções seguintes deste capítulo. 


1.2.3 Sistema em malha aberta 


Feitas as observações iniciais sobre o conceito de sistema, vamos, agora, aprender sistemas em 
malha aberta e em malha fechada. 
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Define-se como sistema em malha aberta o sistema em que não há realimentação da saída 
à entrada, ou seja, a saída é observada apenas a partir de suas amostras instantâneas, não estabe- 
lecendo qualquer conexão com o sinal de entrada, direta ou indiretamente. A Figura 1.2 ilustra grafi- 
camente um sistema em malha aberta. 


Ruído ou 
perturbação 
Sinal de Sinal de 
entrada saída 
Planta 


Figura 1.2 - Sistema em malha aberta. 


Conforme observado na Figura 1.2, um sistema em malha aberta é afetado por perturbações 
cujos efeitos são observados no sinal de saída, não havendo como reduzirem-se tais efeitos a partir 
da mera observação dessa saída. Neste caso, a saída fornece ao observador apenas uma indicação da 
resposta do sistema, não viabilizando meios para o controle automático dessa resposta. 


1.2.4 Sistema em malha fechada 


Define-se como sistema em malha fechada aquele em que há realimentação da saída à entrada, 
ou seja, a saída não somente é observada apenas a partir de suas amostras instantâneas como 
também exerce efeito sobre a entrada do sistema, por meio da geração de um sinal de erro. A Figura 
1.3 ilustra graficamente um sistema em malha fechada. 


Ruído ou 
perturbação 
Sinal de Sinal de 
entrada saída 
Planta 


Sinal de 
realimentação 


Unidade de medição 


Figura 1.3 - Sistema em malha fechada. 


Conforme observado na Figura 1.3, um sistema em malha fechada, ao ser afetado por pertur- 
bações cujos efeitos são observados no sinal de saída, admite o ajuste do sinal de entrada a partir da 
geração de um sinal realimentado que é adicionado (realimentação positiva) ou subtraído (realimen- 
tação negativa). Deste modo, um sinal de erro, obtido da conjunção dos sinais de entrada e de reali- 
mentação, é reinserido continuamente na planta, de modo que a comparação contínua entre o valor 
atual de saída e o valor de setpoint conduz a uma ação de correção que produz uma convergência 
de um valor a outro. Nessas condições, o elemento que gera o sinal de erro a ser inserido na planta 
recursivamente recebe também a denominação de elemento final de controle. 
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Um sistema de controle em malha fechada pode ser mais bem caracterizado pelo desmembra- 
mento da malha em mais subpartes, como ilustrado na Figura 1.4. Nesta concepção, são definidas 
unidades de comparação, de controle em malha direta, de atuação e de medição (ou sensor). Admite-se 
também que haja ruído no processo de medição, igualmente afetando o desempenho da planta. 


Unidade de Ruído ou 
comparação perturbação 


Sinal de 
saída 


Sinal de 
entrada 
Unidade de Unidade de 


E Planta 
controle atuação 


Sinal de 
realimentação 


Unidade de medição 


Ruído de 
medição 


Figura 1.4 - Sistema em malha fechada - descrição detalhada da malha de controle. 


A Tabela 1.1 sumariza vantagens e desvantagens dos sistemas em malha aberta e malha fechada. 


Fique de olho! 


Em rápida análise, sistemas de controle em malha fechada tiram proveito do fato de que a realimentação torna a resposta 
do sistema relativamente insensível a perturbações externas e variações paramétricas do sistema, a um menor custo. Isso 
não ocorre com sistemas de controle em malha aberta, que, por não considerarem a realimentação como proposta, acar- 
retam ao sinal de saída do sistema uma total vulnerabilidade a perturbações externas - sendo assim destinados somente 
a Situações em que os sinais de entrada são previamente conhecidos e imunes a distúrbios externos. 


Tabela 1.1 - Vantagens e desvantagens das categorias de sistema de controle - malha aberta versus malha fechada 


Categoria de sistema Vantagens Desvantagens 


São de simples construção e de fácil 
manutenção. 


» Distúrbios e mudanças na calibração 
provocam erros, e a saída pode apre- 
sentar distorções em relação ao padrão 
esperado. 


São menos dispendiosos que um sistema 
correspondente em malha fechada. 


Malha aberta Não apresentam problemas de estabi- 


est: Para que a saída mantenha a qualidade 


requerida, é necessária uma regulagem 
periódica. 


Sáo apropriados quando há dificuldades 
de avaliacáo dos sinais de saída ou 
quando a medicáo precisa da saída náo 
é economicamente possível. 


continua 
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Categoria de sistema Vantagens Desvantagens 


» O emprego de realimentação faz com 

que a resposta do sistema seja relativa- 

mente insensível a distúrbios e variações 

internas nos parâmetros do sistema. » São de mais complexa implementação e 
de maior custo. 


É possível a utilização de componentes 


Malha fechada relativamente imprecisos e de baixo 
custo para se atingir o controle preciso 
de determinado sistema. 


>» Atrasos temporais podem ocasionar 
oscilação na saída e até mesmo instabi- 
lidade. 

São mais vantajosos nos casos em que 

houver distúrbios e/ou alterações não 

previsíveis nos componentes do sistema. 


Fique de olho! 


A descrição e exemplificação de sistemas de controle sob tal enfoque será revisitada no Capítulo 3, que descreve em 
domínios matemáticos (do tempo e da variável transformada) sinais e sistemas, e mais bem detalhada no Capítulo 5, 
quando se descrevem diferentes estratégias de controle. 


Exercício resolvido 


Considere o sistema de controle de velocidade de um automóvel ilustrado na Figura 1.5. Iden- 
tifique os elementos da malha de controle atuante nesse sistema. 


Solução 


Trata-se de um sistema de controle em malha fechada. A planta é compreendida pelo automóvel 
em si. O sinal de entrada ou de referência é a velocidade desejada. O sinal de saída ou observado é a 
velocidade atual. A unidade de medição do sensor é o velocímetro. O sinal realimentado (ou medi- 
do) é a velocidade medida. O sinal de controle derivado do controlador é a inclinação do pedal do 
acelerador. A unidade de atuação é o motor do veículo. O ruído ou perturbação atuante sobre a 
planta é a inclinação da estrada, e o ruído de medição é a trepidação do ponteiro do velocímetro. 


Inclinação do pedal Inclinação da 
do acelerador estrada 
Velocidade Velocidade 
desejada atual 


Controlador Automóvel 


Velocímetro 


Velocidade 

medida 
Trepidação 
do ponteiro 


Figura 1.5 - Sistema de controle de velocidade de um automóvel. 
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Em suma, como será detalhado em capítulos posteriores e como decorrente da caracterização 
aqui apresentada acerca de sistemas em malha fechada (ou realimentados), um sistema dessa catego- 
ria deve atender a requisitos específicos para sua ação adequada: 


» Robustez a perturbações ou distúrbios: um sistema em malha fechada deve prever e combater 
distúrbios atuantes sobre a planta, de modo a combater com êxito o efeito de tais distúrbios. 


» Erro em regime estacionário: um sistema em malha fechada deve manter controlado o 
erro quando da convergência do sinal de saída ao valor desejado. 


» Características em regime transitório: um sistema em malha fechada deve atender a 
requisitos específicos em sua saída enquanto esta não atinge o valor de setpoint, condição 
que recebe a denominação de regime transitório. Valores específicos e desejados de tempo, 
bem como valores máximos que a saída deve apresentar nessa condição de regime, são 
parâmetros de projeto de controladores, como será visto nos Capítulos 4 e 5. 


1.3 Modelagem de sistemas físicos 


Na descrição de sistemas, o número elevado de variáveis envolvidas torna complexa sua pre- 
cisa caracterização, conquanto as respostas a determinados sinais de entrada possuam dependência 
com uma enorme gama de parâmetros e critérios. Nesse sentido, a tarefa de modelar um sistema, 
restringindo suas características a um menor, mas suficiente conjunto de descrições, torna sua repre- 
sentação matemática adequada ao estudo, mantendo ainda a robustez necessária para predizer seu 
comportamento a praticamente todas as possíveis condições de entrada. 


Ao se definir assim um modelo como uma representação substitutiva da realidade de um sis- 
tema, estabelece-se uma formulação ou equacionamento que expressa as características de um sistema 
físico, em termos matemáticos. Salienta-se, porém, que, por mais precisão que um modelo físico 
contenha, a partir de equações cada vez mais complexas, esse modelo nunca será exato. 


Fique de olho! 


Um modelo corresponde à descrição matemática das características dinâmicas de um sistema adequado à sua análise. A 
partir de seu modelo, um sistema pode ser avaliado por meio de ferramentas computacionais, advindo daí a necessidade 
de que tal modelo seja simples o suficiente para ser implementado, mas complexo o bastante para garantir precisão em 
seus resultados. 


Sistemas podem ser divididos, quanto ao número de sinais de entrada e correspondentes sinais 
de saída, em duas categorias: SISO (Single Input - Single Output), Figura 1.6, ou monovariáveis, e 
MIMO (Multiple Input - Multiple Output), Figura 1.7, ou multivariáveis, de acordo com suas caracte- 
rísticas de mapeamento de entradas e saídas. Podem ainda ser contínuos, se existem para todo valor 
real de tempo, ou discretos, se operam apenas sobre valores inteiros de tempo. 


x(t) y(t) 
— Sistema ——> 


Figura 1.6 - Sistema SISO. 
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Xx (t) 


Sáo variados os exemplos de modelos em engenharia, como mostram as Figuras 1.8 a 1.10. 


Sistema 


HA nl) 
>> y alt) 


m alt) 


HH yli) 


Figura 1.7 - Sistema MIMO. 


= 


Figura 1.8 - Modelo de uma fonte de tensão real E. 


Figura 1.10 - Modelo no domínio transformado de sistema realimentado para controle de nível em um tanque. 


Para a obtenção de modelos de sistemas, são utilizadas leis físicas que os descrevam, usual- 
mente baseadas em equações diferenciais lineares, a partir de que saída é obtida para diferentes 
possibilidades de entrada: resposta a um sinal impulsivo (no domínio do tempo), resposta a uma 
entrada senoidal (no domínio da frequência), entre outras. 
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k e 
F(t) 


Figura 1.9 - Modelo de um sistema massa-mola-amortecedor. 


Desta forma, o uso conjugado de adequadas técnicas que descrevam o comportamento dos sis- 
temas a entradas específicas em regime transitório e permanente, aliado a um bem descrito modelo 
matemático que traduza com propriedade as características físicas reais de uma planta, facilita o 
estudo de estratégias de controle aplicadas a tais plantas. Os próximos capítulos descreverão com a 
devida base os procedimentos dedicados ao projeto de sistemas de controle, partindo-se da neces- 
sária base matemática à quantificação específica de parâmetros de controladores. 


Vamos recapitular? 


Foram descritas neste capítulo as noções básicas de controle de sistemas/plantas representativos 
de fenômenos naturais, bem como ressaltada a importância do modelamento de tais sistemas, visando 
ao adequado projeto de controladores que se prestem a ajustar sinais de saída desses sistemas a valores 
desejados. 


O capítulo abordou, ainda, características detalhadas de sistemas de controle em malha aberta e 
malha fechada, de modo a pontuar vantagens e desvantagens de tais estratégias, úteis às etapas de projeto 
de controladores que são descritas nos próximos capítulos. 


O Agora é com você! 


1) Destaque e caracterize os elementos de controle envolvidos em um processo simples 
de troca de calor existente em um freezer doméstico. Aponte eventuais modalidades de 
ruído presentes na planta e na medição de temperatura, apontando a função de cada 
elemento de controle. 


: 2) Identifique nos exemplos a seguir a aplicação de sistemas de controle de malha aberta 

: ou malha fechada: (a) controle de uma máquina de secar roupas; (b) controle manual 
de abertura de uma válvula de chuveiro, de acordo com sua temperatura; (c) controle de 
trajetória de um foguete. 


: 3) Para você, o corpo humano é um sistema de controle em malha aberta ou malha 
fechada? Justifique sua resposta a partir de exemplos. 
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Elementos de 
Algebra Linear 


Para começar 


Este capítulo tem por objetivo apresentar o ferramental matemático mínimo necessário ao enten- 
dimento de como sistemas de controle são aplicados a outros sistemas, com o objetivo de adequar seus 
sinais de saída a requisitos desejados. A álgebra linear se refere ao campo da matemática que estuda de 
modo detalhado os sistemas lineares, compostos de equações algébricas ou diferenciais. Como a tais sis- 
temas são aplicados sistemas de controle para obtenção de sinais de resposta específicos, o entendimento 
das relações matemáticas em álgebra linear é relevante ao bom projetista de controle. 


Serão apresentadas neste capítulo entidades matemáticas como vetores, matrizes e sistemas, tra- 
zendo o embasamento necessário aos capítulos posteriores. 


2.1 Vetores e matrizes 


Define-se vetor como o segmento de reta orientado dotado de módulo (ou magnitude), dire- 
ção e sentido, representando um espaço vetorial de n dimensões. Aos vetores podem ser aplicadas 
operações de adição, subtração, escalonamento e inversão. A eles também se aplicam propriedades 
de comutatividade, associatividade e distributividade. 


Um vetor vV é representado a partir de suas componentes em n dimensões. No espaço tri- 
dimensional, obedece-se a notação de decomposição de um vetor em três dimensões, em que 
se dispõem três vetores unitários nas direções X, Y e Z, denominados í SÉ). Por exemplo, o vetor 
y= di +3]+6k é descrito no espaço vetorial tridimensional, possuindo componentes escalares 


4,3 e 6, ao longo dos eixos X, Y e Z, respectivamente. Deste modo, tais componentes descrevem as 
coordenadas cartesianas tridimensionais da extremidade do vetor em questão, como na Figura 2.1. 


A A A 
— (4i + 3j + 6k) 
E O 
Figura 2.1 - Vetor V = 4i +3j+ 6k . 


A definição de matriz é uma extensão da definição da entidade vetor, representando uma 
tabela de m linhas e n colunas (e assim sendo rotulada de matriz de dimensões mxn ), também útil 
na análise de sistemas de equações lineares. Quando uma das dimensões da matriz é igual à unidade, 
ela se reduz a um vetor. São exemplos de matrizes: 


1 
2 3 
M= 12 , A =| 4 | , entre outros. 
2 


Nos exemplos citados, as dimensões são 2x2 e 3x1, respectivamente. Genericamente, o ele- 
mento m; ; de uma matriz se situa na i-ésima linha e na j-ésima coluna. 
2.1.1 Operações com vetores 


Vetores podem ser operados matematicamente sob diversas situações: 
a) Adição/subtração de vetores 


A adição/subtração de vetores obedece à adição/subtração de seus componentes ao 


longo de cada eixo cartesiano. 


a AA A 
Ted + E] 


= A A A 
W=w,1 + w,J+w,k 
DEW =(v, +w, i+ (v, +w,)J+(v, + w,)k (2.1) 


E S A A A 
V-W=(v, =, =w +, —W,)k 
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b) Multiplicação de um vetor por um escalar 


A multiplicação de um vetor por um escalar é similar à operação em escalares puros. 


= AAA 
v= itv J+r E 


W=avV , AER (2:29) 
z A A 
W=av,i+0v,j+0v,k 


> 


c) Produto escalar ou produto interno 


O produto escalar (ou interno) entre dois vetores é expresso pela Equacáo 2.3, resul- 
tando em um escalar. 


z ^ A A 
veniti Jen 


x A A A 
= mem] (2.3) 


VW = VW F VyW HNW, 


y 
Valem para o produto escalar as propriedades dispostas na Tabela 2.1. 


Tabela 2.1 - Propriedades do produto escalar entre dois vetores 


Propriedade Operação 


Elemento nulo 


d) Produto vetorial 


O produto vetorial entre dois vetores é expresso pela Equação 2.4, resultando em um 


vetor (admite-se que tais vetores estejam definidos no espaço tridimensional R?). 


m A AA 
fa a ] nd 


A A 
j k (2.4) 
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Na equação anterior, |e| simboliza a operação determinante de uma matriz, gerando como 
resultado da operação de produto vetorial um vetor perpendicular aos vetores operados. Valem para 
o produto vetorial as propriedades dispostas na Tabela 2.2. 


Tabela 2.2 - Propriedades do produto vetorial entre dois vetores 


Propriedade Operação 


XxW=-WxV 


w) 


Conversão em vetores com produtos escalares | (Vxw)xU =(V.U).w —(V.w).ú 


O produto escalar e o módulo do produto vetorial entre dois vetores podem ainda ser expres- 
sos em função do ángulo O formado entre tais vetores, como mostram as Equações 2.5 e 2.6. 


S ES 
VETA, IA 
= A A A 
w=w,i+w,j+w,k (2.5) 


v.w =|v||w|cos 0 

= A A A 

Fendi, rã 
A 


w=w i+w J+ w,k (2.6) 


|V x w| = |v||w| seno 


Uma operação interessante ainda no campo matemático dos vetores é o chamado produto misto, 
expresso pela Equação 2.7, que relaciona três vetores e uma diferente combinação de produtos entre 
eles. Tal operação é bastante utilizada quando da decomposição de vetores em distintas dimensões. 


= A AA 
= TE] EE 


A A A A 
w=w,i+w, j+w,k 
E A AA 
u=uçi+u,j+u,k 
27, 
la U, U, (2.7) 


(Vxw)ui=|v, Ye Ve 


Wi Wy Wa 
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2.1.2 Operações com matrizes 


Matrizes podem ser operadas matematicamente sob diversas situações: 


a) Adição/subtração de matrizes 


A adição/subtração de matrizes obedece à adição/subtração de seus componentes de 


mesmo índice. 


a1 

A =| az 
a31 

by 
B=|b,; 
bs; 

aj +b; 
A+B=|a,,+b,; 
as, +Dy, 
a Dj, 
azı ~ b31 


a 3 

do “3 

a32 3 

bo bj 

by by 

by bas | 

ap +b? ag+by (2.8) 
az +b az +bz3 
a32 +by asa + b33 
aby, abs 
a,b, az -b23 
a3, bz a33- b33 


b) Multiplicação de uma matriz por um escalar 


A multiplicação de uma matriz por um escalar é similar à operação em escalares puros. 


ajy Mg 
A=|a a, 
a3; d32 
B=aA 

da, da 
B=|aa,, cas 
Os, E) 


c) Multiplicação entre duas matrizes 


a13 
a23 
a33 


(2.9) 
da, 
da, 


Oass 


A multiplicação entre duas matrizes é bem definida sob a condição de que o número 
de colunas da primeira matriz seja igual ao número de linhas da segunda matriz, 
ou seja, para que seja válido o produto entre duas matrizes A e B, é necessário 
que possuam respectivamente dimensões mxn e nxp, tendo a matriz resultante 
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do produto dimensões mxp. Sob tais condições, cada elemento da matriz produto 
C= AB é dado por: 


Ee aiibi j +a;ob, j +..+a; b (2.10) 


J in n,) 


(9) Exercício resolvido 


1-2 
1) Determine o produto entre as matrizes A = E | e B= 


Solucáo 


A matriz produto possui 2 linhas e uma coluna: 


mi2 
AB= 
[3 4 a 
[14+1.0+2:(=1) 
AB= 
[3.44+4.0+1(-1) 
le 
AB= 
al 


O resultado obedece, assim, às dimensões esperadas para o resultado da operação produto 
entre matrizes. 


2.1.3 Matrizes especiais 


Algumas matrizes recebem denominações especiais, visando a facilitar algumas operações 
nesse conjunto de entidades matemáticas. 


a) Matriz identidade 


A matriz identidade 1 é a matriz quadrada de dimensões n x n em que a diagonal 
principal é formada por algarismos 1 e os demais elementos são nulos. 
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[100 0 

010 0 

0 0 se 0 
pE ......0 (2.11) 

0 

s mm E 

[0000001 


b) Matriz transposta 


Define-se matriz transposta de uma matriz conhecida M àquela obtida da inversão 
linha/coluna de seus elementos m; j> tal que: 


i us 
m ja Mji 


Importante observar que (A BC.D....Z)T = ZT....CT.BT.AT. 


Exercício resolvido 
tod 2 
2) Determine a matriz transposta de A = p | 


Solução 


c) Matriz simétrica 


Define-se matriz simétrica de uma matriz M aquela que é igual à transposta de M. Tem- 
-se como condição necessária para tal que M seja quadrada. A matriz simétrica é útil na 
determinação reversa de respostas a partir de matrizes representativas de sistemas. 


d) Matriz inversa 


Define-se matriz inversa M~! de uma matriz quadrada M aquela tal que 
MM! =M~!.M =1. A condição necessária e suficiente para que uma matriz qua- 
drada M seja inversível é que det(M) * 0. 
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(O) Exercício resolvido 
Lo o 


3) Verifique sea matriz A=|2 —1 0| possui inversa. 


4 3 1 
Solução 
o a 
de(A)=2 -1 0|=-1+12+0-(-8+4+0)=11-(—4)=15+0 
4 3 1 


A Tabela 2.3 sumariza propriedades das matrizes quanto á sua classificacáo. 
Tabela 2.3 - Propriedades das matrizes quanto à sua classificação 


Tipo de matriz É quadrada? É inversível? Qual sua transposta? 


Matriz identidade I,, Sim, ela própria Ela própria 


La Es e 


2.2 Autovalores e autovetores 


Considere-se uma transformação linear de um vetor em outro, T: V > W tal que 


T(v)=Av, (2.12) 


em que À é um escalar dito autovalor e v é dito autovetor, se v > 0. Uma vez que uma transforma- 
ção linear pode ser escrita pelo produto de uma matriz por um vetor, pode-se escrever: 


T(v)=Av. (2.13) 
Estabelecendo-se igualdade entre as equações anteriores, pode-se escrever: 
(A-ADv=0, (2.14) 


onde A é uma matriz quadrada. Os vetores v > 0 para os quais há um valor de À que seja solução da 
Equação 2.14 são denominados autovetores da matriz 4, e À são denominados autovalores associa- 
dos àqueles autovetores. Para que isso seja possível, é condição que det(A —AI) = 0, o que gera um 
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polinômio de n-ésimo grau denominado polinômio característico (a ser abordado no próximo capí- 
tulo), cujas raízes são exatamente os autovalores da matriz A. 


(9) Exercício resolvido 
à transformação linear T:R? —> R? , expressa por 


4) Obtenha os autovalores associados 
(x,y) (4x +5y,2x + y). 


Solucáo 


A obtenção de autovalores segue o exposto na teoria supracitada: 


E | f i o 5 | 
det(A-A)=0... EN = =0 
Doi 01 2 dd 


(4-AX1-1)-10=0 


( Vamos recapitular? 


Foram apresentados neste capítulo conceitos fundamentais em álgebra linear, considerados rele- 
vantes às técnicas de controle a serem descritas nos próximos capítulos. 


O capítulo abordou, entre outros elementos, os conceitos e propriedades de vetores e matrizes, 
enunciando também a definição de autovetores e autovalores, úteis em determinadas técnicas de ajuste 


de sinais de saída provenientes de sistemas lineares. 
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O Agora é com você! 


: = A Å A 
: 1) Obtenha os produtos vetorial e escalar entre os vetores V=31+4j+5k e 


A A A 
w=21+2)+4k. 
3 2 3 
: 2) Quala matriz inversa da matriz A=|2 1 2|? 
: 141 


O que se pode observar a respeito de tal matriz? 
: 3) Obtenha os autovalores associados à transformação linear T : R? — R2, expressa por 
(x, y) > Gx T y,2x E y). 


: 4) Pesquise sobre matrizes adjuntas e obtenha sua relação com o cálculo da matriz 
inversa de uma matriz qualquer. 
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Sinais e 
Sistemas Lineares 


Para começar 


Este capítulo tem por objetivo apresentar os conceitos de sinais e sistemas lineares, dentro do 
escopo da Teoria de Controle. O capítulo traz representações de sinais e sistemas em domínios distintos, 
abordando ainda a caracterização de sistemas em termos de raízes de suas funções de transferência e de 
sua estabilidade em regime permanente. 


As informações contidas neste capítulo abrem as portas para a escolha adequada de que modali- 
dade de projeto de sistema de controle melhor se encaixa como solução. 


3.1 Conceitos básicos e propriedades 


Como visto no Capítulo 1, um sinal é uma representação matemática de um fenômeno da 
natureza. Sob o ponto de vista da Teoria de Controle, são relevantes os sinais que alimentam siste- 
mas lineares, no sentido de que a matemática envolvida é facilmente mapeada e desenvolvida por 
meio de ferramentas de cálculo diferencial e de álgebra linear, além de se dispor de um sem número 
de técnicas de matemática aplicada - caso das transformações de domínio -, assegurando uma 
caracterização completa de sinais e sistemas em termos de equações de fácil resolução. 


Nesse sentido, alguns conceitos podem ser estabelecidos, a partir do conceito primário de sistema. 
a) Sistema causal (ou não antecipativo) 


É todo sistema que responde somente a partir do instante em que é excitado por uma 
entrada, ou seja, depende de amostras de entrada atuais ou passadas, somente. Se y(t, T) 


é a saída no instante t de um sistema excitado por um sinal de entrada aplicado no 
instante T, isto significa que y(t,t)=0 para t< t. 

Exercício resolvido 

1) Verifique se o sistema descrito pela equação y(t) = x(t) +2x(t+2) é causal. 

Solução 

Como a saída y(t) é função de componentes de entrada que ocorrem em instantes de tempo 


futuros, a saída do sistema se antecipa à entrada, de modo que o sistema é não causal. 


b) Sistema invariante no tempo 


É todo sistema que possui parâmetros constantes, dependendo somente das variá- 
veis de entrada, e jamais do tempo. Em outras palavras, se um sistema possui resposta 
y(t) a um sinal de entrada x(t), a aplicação de um sinal x(t-t) conduzirá a uma 
saída y(t—7). 


Exemplo: 
y(t) = x(t) + x(t-2) > é invariante no tempo. 


y(t) = x(3t) + x(t—1) > é variante no tempo. 


3.2 Representacáo no domínio do tempo 


A representação de sinais e sistemas no domínio do tempo permite a análise sob o ponto de vista 
da operação de convolução entre sinais, para se calcular a saída de um sistema a uma dada entrada. 
Seja x(t) um sinal arbitrário aplicado à entrada de um sistema linear e h(t) a resposta de tal sistema a 
uma entrada impulsiva (t). O sinal de saída y(t) para tais condições pode ser expresso por: 


oo 


y(0)=x(0+h(0)= f x(1)h(t-0)d (3.1) 


—00 


A integral expressa na Equação 3.1 é dita integral de convolução, e por vezes é de difícil solução 
analítica, recorrendo-se assim (a) à solução gráfica (que foge do escopo deste livro), Figura 3.1, ou 
(b) à representação de sinais em um domínio alternativo transformado, como o da transformada de 
Laplace, expresso na próxima seção deste capítulo. 


x(t) y(t 
h(t) 


Figura 3.1 - Sistema linear caracterizado por sua resposta ao impulso h(t). 
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3.3 Representação no domínio da frequência - 
a transformada de Laplace 


A representação de sinais e sistemas lineares no domínio da transformada de Laplace assegura 
um modo simples e compacto de se avaliar tais entidades, possibilitando ainda que se possam ana- 
lisar tais sistemas sob o enfoque de equações algébricas em vez de equações diferenciais, usuais na 
representação temporal. 


De modo generalizado, a transformada de Laplace admite ainda o emprego de tabelas de pares 
de transformação direta e inversa e de propriedades (constantes dos Apêndices C e D desta obra), 
conduzindo a uma análise elegante e precisa, além de possuir uma forte associação da caracterização 
das respostas à posição das raízes das equações algébricas associadas aos sistemas. 


A transformada de Laplace X(s) de um sinal x(t) é dada por 


00 


X(s) = feat (3.2) 
0 


Nesse sentido, sinais e equações de elaborada caracterização são facilmente transportados, no 
domínio da transformada de Laplace, para equações polinomiais em s. Merece destaque a defini- 
ção da variável complexa s = o + jæ, que é nada mais que um caso generalizado da transformada de 
Fourier. O estudo da transformada de Fourier não se situa no escopo deste livro, mas é recomendado 
igualmente na análise de sinais, podendo constituir estudo complementar ao aluno. 


(9) Exercício resolvido 


y 120 
2) Determine a transformada de Laplace do sinal degrau unitário x(t) = u(t) = E o 
t< 


Solução 
O cálculo da transformada de Laplace obedece à fórmula de transformação, considerando que 
o sinal analisado, degrau unitário, não existe para instantes negativos de tempo. Assim: 


X(s)= f x(De-stdt 
0 


X(s) = f Le-Stdt 
0 


00 


X(s) = = em 
s 


0 


Rea 
S 
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3.4 Caracterização de polos e zeros 


Um modo imediato de se avaliar o comportamento de um sistema a uma dada entrada con- 
siste em se analisar as raízes da transformada de Laplace de sua resposta a uma entrada impulsiva. 
Seja H(s) a chamada função de transferência de um sistema, equivalente à transformada de Laplace 
da resposta ao impulso h(t) desse sistema, e suponha-se que H(s) é racional, na forma 


N(s) 
H(s)=—— > 

D(s) 
em que as raízes dos polinômios do numerador e denominador são definidas, respec- 
tivamente, como zeros e polos. A localização das raízes do termo D(s) , ou seja, os polos do sistema, 
fala muito sobre o comportamento desse sistema quando sujeito a uma entrada arbitrária. 


Genericamente, os polos de H(s) são complexos, e a análise sobre a existência de partes real e 
imaginária desses polos permite que se conclua sobre a natureza da resposta do sistema. Antes de se 
tecerem conclusões acerca de tais características, o exemplo a seguir de cálculo de transformada de 
Laplace é bastante interessante para a análise aqui realizada, sendo fundamental ao entendimento 
apurado de controladores em outro domínio diferente do temporal. 


(9) Exercício resolvido 
e t>0 


3) Determine a transformada de Laplace do sinal x(t) = e-“tu(t) = N 
tg 


Solução 


Tal como no exercício anterior, a transformada de Laplace do sinal analisado é dada por: 


X(s) = Í x(t)e-stdt 
0 
X(s) = je aL 
0 
X(s)=—est| = (e —e0) 
s+0 | sta 
o 
s+a 


Como observado no exemplo anterior, existe uma relação direta entre a posição de um polo 
a e sua resposta temporal correspondente e-ºt. Desta maneira, fica nítido o comportamento de tal 
resposta, se convergente ou divergente, a partir da natureza do polo a. No caso de a complexo, a 
exponencial temporal se tornaria também complexa, o que induz a um comportamento senoidal, 
portanto oscilatório, da resposta. 
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A Tabela 3.1 mostra as características de resposta de um sistema a uma entrada em degrau 
unitário, de acordo com a característica dos polos de tal sistema (para efeito de polos variados, 
toma-se como predominante o polo mais à direita no plano complexo). 


Tabela 3.1 - Propriedades da resposta de um sistema de acordo com a localização de seus polos 


Característica 
da resposta 


RR o mea pe i Ai 


Parte real Parte imaginária Exemplo 


aa < i cli 


3.5 Estabilidade de sistemas lineares 


les ji dl is 


A condição de estabilidade de um sistema linear é provavelmente a mais importante a ser 
assegurada no comportamento de um sistema. Um sistema admite em sua resposta, sob certas con- 
dições, até mesmo um comportamento oscilatório, desde que convirja a determinado valor após 
algum intervalo aceitável de tempo. Todavia, um sistema monotonicamente instável jamais será 
desejável, conquanto haja especificações de regime permanente que devem ser atendidas. Essa talvez 
seja a maior finalidade de sistemas de controle: prover a sistemas teoricamente instáveis condições 
para que atinjam a condição de estabilidade. 


O próximo capítulo versará sobre a estabilidade de sistemas lineares, especificamente para 
cada categoria de sistema estudado. 


Amplie seus conhecimentos 


Um exemplo clássico de sistema instável é o apre- 
sentado nas figuras a seguir. A primeira ponte sobre o 
estreito de Tacoma, em Puget Sound, Washington, 
mostrada na Figura 3.2, foi aberta ao tráfego em 1º de 
junho de 1940. Após algum tempo, percebeu-se que a 
ponte oscilava após ventanias. Cerca de quatro meses 
após sua inauguração, uma ventania produziu uma 
oscilação que cresceu em amplitude até a ruptura total 
da ponte, que pode ser vista na Figura 3.3. Este é um 
exemplo típico de sistema instável que, após a aplicação 
de uma entrada específica, se põe a oscilar até atingir 
valores divergentes do equilíbrio. 


Autor desconhecido/Wikimedia Commons 


às 


Figura 3.2 - Tacoma Narrows Bridge antes da ruptura. 
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Figura 3.3 - Tacoma Narrows Bridge após a ruptura. 


_ Vamos recapitular? 


Foram apresentados neste capítulo os conceitos de sinais e sistemas lineares, do ponto de vista de 
suas possíveis representações nos domínios do tempo e da transformada de Laplace, que é de caráter fre- 
quencial. Por meio do entendimento das raízes do sistema no domínio da frequência, é possível estimar 
o comportamento dos sinais de saída desses sistemas, quando sujeitos a sinais de entrada elementares. 


O capítulo abordou, entre outros elementos, os conceitos de estabilidade e oscilação de sinais de 
saída, elencando-os à localização de polos e zeros no plano complexo. 


: Agora é com você! 
1) Utilizando os pares de transformada de Laplace (Apêndice), determine a transforma- 
da de Laplace dos sinais x(t) = e3tu(t), x(t) = cos(2t)u(t) ex(t) = třu(t). 


2s+1 


s? +6s+8 
Tal sistema é estável ou instável? Oscilatório ou não? Justifique sua resposta. 


: 2) Identifique os zeros e polos da função de transferência H(s) = 


: 3) Explique com suas palavras o fenómeno da instabilidade observado em um carro que 
desce sem freios uma ladeira com 45º de inclinação. 
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Sistemas de 
Primeira e 
Segunda Ordens 


Para começar 


Este capítulo tem por objetivo apresentar sistemas lineares de primeira e segunda ordens, con- 
siderados os mais comuns e relevantes como modelos no estudo de fenômenos realísticos modelados 
matematicamente. 


As informações contidas neste capítulo constituem uma importante base para a compreensão 
aprofundada de técnicas de projeto de sistemas de controle, abordadas a partir do próximo capítulo. 


4.1 Estabilidade de sistemas de primeira 
e segunda ordens 


A estabilidade de sistemas remete a um grande número de problemas em que é preciso man- 
ter constante uma dada variável de saída. Em alguns casos, é desejável um transitório acelerado, ao 
passo que em outros casos é preciso minimizar o erro de regime permanente. Nesses casos, o pro- 
jetista deve saber projetar adequadamente os filtros de controle do sistema. O primeiro passo, no 
entanto, é saber especificar matematicamente os índices de desempenho desejados para a resposta. 


Nesse contexto, a estabilidade de sistemas de primeira ordem é de simples compreensão, pois 
esses sistemas apresentam diversos polos que invariavelmente produzem respostas exponenciais 


oscilatórias ou não. A forma de uma função de transferência composta de termos em primeiro grau 
é dada por: 


(4.1) 


H(s) = 
s+0 


Sob tais condições, a saída será composta de uma função exponencial, convergente ou não, a 
depender da polaridade do sinal da parte real de cada polo. 


Para sistemas de segunda ordem, a análise é mais complexa. Usualmente, um sistema de 
segunda ordem é representado por: 
o? 


H(9)=>—— |-——— 4.2 
(9) s? +2£0,5+0,? po 


A condição de estabilidade de um sistema linear é, provavelmente, a mais importante a ser 
assegurada no comportamento de um sistema. Um sistema admite em sua resposta, sob certas con- 
dições, até mesmo um comportamento oscilatório, desde que convirja a determinado valor após 
algum intervalo aceitável de tempo. 


Todavia, um sistema monotonicamente instável jamais será desejável, conquanto haja especi- 
ficações de regime permanente que devem ser atendidas. Essa talvez seja a maior finalidade de siste- 
mas de controle: prover a sistemas teoricamente instáveis condições para que atinjam a condição de 
estabilidade. 


4.2 Especificações de desempenho 


São relevantes na observação da resposta de um sistema de segunda ordem alguns parâme- 
tros de projeto, a partir do polinômio característico do sistema (denominador da expressão na 
Equação 4.2): 


» Tempo de subida t, (rising time): intervalo de tempo que a resposta leva entre cruzar pela 
primeira vez um determinado limite inferior e a primeira vez que cruza um determinado 
limite superior. Estes limites são definidos em função de uma porcentagem do valor final 
(por exemplo, 10 - 90%, 5 — 95% e 0 — 100%). 


» Tempo de pico t (peak time): tempo necessário para a resposta atingir seu valor máximo, 
ou seja, seu sobressinal. E modelado por: 


T 
t = 


e 0, y1-8? 


»  Sobressinal SS (overshoot): máximo valor de pico da curva de resposta medido a partir do 
valor unitário. É modelado por: 
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» Tempo de acomodação ou tempo de assentamento t : tempo necessário para a curva de 
resposta alcançar e permanecer dentro de uma faixa em torno do valor final VF, de mag- 
nitude especificada por uma porcentagem absoluta daquele valor final (por exemplo, 95% 
ou 98%). O tempo de acomodação é modelado por: 


3 
t, = FA (critério de 95% da resposta final) 
n 


4 
t, = 20 (critério de 98% da resposta final) 
n 


As especificacóes de desempenho de um sistema linear podem ser caracterizadas graficamente, 
como vemos na Figura 4.1: 
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Figura 4.1 - Caracterização das especificações de desempenho de um sistema linear - resposta ao degrau unitário. 


A partir de tais requisitos, o projetista pode analisar qual a melhor estratégia de controle, de 
maneira a privilegiar um ou outro. Será visto nos capítulos posteriores que a escolha por mais veloci- 
dade de resposta ou por menos sobressinal influi diretamente no modelo de controlador selecionado. 


Um problema de controle se concentra, assim, em determinar um modo de se alterar um sis- 
tema físico conhecido, a fim de atender às especificações de desempenho desejadas. Como o sistema 
físico é usualmente imutável, a solução consiste em se projetar controladores que, associados ao sis- 
tema, possibilitem que os requisitos desejados sejam obtidos. A escolha adequada do controlador a 
ser utilizado admite, assim, diferentes possibilidades de solução, e a melhor é aquela mais adequada 
à necessidade de projeto. 
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Exercício Resolvido 


1) Para o sistema dado por H(s) = = , determine: 


s? +4s+13 
seus polos, seu tempo de pico, seu sobressinal percentual e seu tempo de acomodação. 
Solução 


A obtenção das raízes do polinômio característico é obtida por meio da resolução da equação 
de segundo grau associada: 


s =-—2+ j3 


" 


s? +4s+13=0 i 
S=? B 


Do polinômio característico obtém-se œ, = 3,605 rad/s e € = 0,555. 
O tempo de pico é calculado como: 


poe z ==1,047s 


PO il= so 
Ae 


O sobressinal é calculado como: 


fe k 
SE =0,234 


O tempo de acomodação é calculado como: 


4.3 Sistemas de primeira ordem 


A função de transferência de um sistema de primeira ordem, como enunciado na Equação 4.1, 
é dada por 


H(s) = É, (4.3) 
S+ 
o que conduz à resposta ao impulso 
h(t) = ke-tu(t), (4.4) 
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o que requer a >0 para que o sistema seja estável. O termo T=1/a é usualmente denominado 
constante de tempo do sistema. Deste modo, pode-se reescrever: 


k ok 


H(s) = = 
(9) s+1/tT 1s+1 


(4.5) 


que corresponde à equação diferencial 


dy E 4.6 
a +ay(t) = kx(t). (4.6) 


A Figura 4.2 mostra a resposta de um sistema de primeira ordem ao impulso. 


A resposta de um sistema linear invariante no tempo de primeira ordem ao degrau unitário 
u(t) é dada por: 


s(t) = (1-e-“Du(t). (4.7) 


0 T 21 31 t 


Figura 4.2 - Resposta ao impulso de um sistema de primeira ordem. 


A Figura 4.3 mostra a resposta de um sistema de primeira ordem ao degrau. 


T t 


Figura 4.3 - Resposta ao degrau de um sistema de primeira ordem. 


Há de se avaliar como parâmetros relevantes na resposta de sistemas de primeira ordem os seguintes: 


» Tempo de subida: 
t 10-90%) = 2,27 
t (595%) = 2,957 
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» Tempo de estabilização: 
t (5%) = 31 
Quanto ao tempo de subida, é infinito, já que a resposta nunca ultrapassa o valor final antes do 


infinito. Uma maior constante de tempo conduzirá a um maior intervalo de tempo até que o sistema 
se estabilize, condigáo em que o polo se situa cada vez mais próximo do eixo imaginário. 


4.4 Sistemas de segunda ordem 


Como visto anteriormente, um sistema de segunda ordem é representado por: 


2 
O 


H(s) = (4.8) 


s? +260,5+0,? * 


em que £ é o chamado fator de amortecimento do sistema e (>, é a chamada frequência natural 
do sistema. 


A Figura 4.4 mostra a localização de polos complexos conjugados de um sistema de segunda 
ordem, no plano complexo. Os polos indicados na figura são -Cw, + jw,1-C?, ou, em forma 
compacta, —1/1, tjw,, onde 1, =-1/50, e q, =0,y1-8?. 


jon 


On 


Figura 4.4 - Resposta ao degrau de um sistema de primeira ordem. 
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Pode-se classificar um sistema de acordo com o valor de seu fator de amortecimento, como 
visto nas subseções a seguir. 
4.4.1 Sistemas subamortecidos: O<f[< 1 


Sistemas de segunda ordem subamortecidos, quando sujeitos a entradas impulsivas, apresen- 
tam como resposta: 


h(t) = (Ae neiet + Aee Int) (4.9) 


A Figura 4.5 mostra diferentes respostas de tais sistemas ao impulso, em função de diferentes 
valores de fator de amortecimento, para uma mesma frequência natural. Considera-se na ilustração, 
além do caso de subamortecimento, o caso crítico de amortecimento (6 = 1). 


07 


0.8- 
0.6- 
0.47 f 


0.24F 


t 
0 2 4 6 8 10 12 14 16 18 20 


Figura 4.5 - Resposta ao impulso de um sistema subamortecido de segunda ordem. Os valores de fator de 
amortecimento estão indicados nas curvas, considerando uma mesma frequência natural. 


Uma análise rápida de tal resposta mostra que oscilações somente acontecem em decorrên- 
cia do segundo termo exponencial, ao passo que o primeiro termo exponencial induz a resposta ao 
amortecimento da amplitude da resposta. 


Os mesmos sistemas, sujeitos a entradas do tipo degrau unitário, apresentam como resposta 
(obtida da integração de h(t)): 


s(t) =|1- La ett, sen(o,t+0) |u(t) (4.10) 


A Figura 4.6 mostra diferentes respostas de tais sistemas ao degrau, em funcáo de diferentes 
valores de fator de amortecimento, para uma mesma frequência natural. Considera-se na ilustração, 
além do caso de subamortecimento, o caso crítico de amortecimento (6 = 1). 
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Figura 4.6 - Resposta ao degrau de um sistema subamortecido de segunda ordem. Os valores de fator de 
amortecimento estão indicados nas curvas, considerando uma mesma frequência natural. 
4.4.2 Sistema sobreamortecido: [> 1 


Sistemas de segunda ordem sobreamortecidos, quando sujeitos a entradas impulsivas, apresen- 
tam como resposta: 


h(t) = M(ebi! —eP2bu(t), (4.11) 
0) 


com M = —=., 
lt 


A Figura 4.7 mostra diferentes respostas de tais sistemas ao degrau em função de diferentes 
valores de fator de amortecimento, para œ, =1 rad/s. 


t 
0 2 4 6 8 10 12 14 16 18 20 


Figura 4.7 - Resposta ao impulso de um sistema sobreamortecido de segunda ordem. Os valores de fator de 
amortecimento estáo indicados nas curvas, considerando frequéncia natural unitária. 


Sistemas de Controle - Matemática Aplicada a Projetos 


Os mesmos sistemas, sujeitos a entradas do tipo degrau unitário, apresentam como resposta 
(também obtida da integração de h(t)): 


s(t) = a T u(t) (4.12) 
Pı P2 


A Figura 4.8 ilustra o comportamento desta resposta para frequência natural unitária, para dis- 
tintos valores de fator de amortecimento, incluindo o caso de amortecimento crítico. Como obser- 
vado, a resposta é sempre crescente, tendendo assintoticamente para o setpoint, sem oscilações. Tal 
como no caso de subamortecimento, a frequência natural aparece somente como um fator de escala 
no domínio do tempo (eixo horizontal). Fica também evidente que a primeira exponencial é domi- 
nante em relação à segunda, aproximando-se do caso de um sistema de primeira ordem com polo p}, 
de modo mais nítido na resposta ao degrau do que na resposta ao impulso, pelo fato de que ocorre 
um escalonamento da amplitude da segunda exponencial por p,, maior em módulo do que p,. Obvia- 
mente, o sobressinal da resposta ao degrau para sistemas sobreamortecidos é nulo, pois a resposta 
nunca ultrapassa o valor máximo. Já aproximações para os tempos de subida e de estabilização 
podem ser facilmente obtidas na condição de que £ >> 1, mantêm-se as expressões: 


» Tempo de subida 
t,(10-90%) = 2,21 
t (5-95%) = 2,951 


» Tempo de estabilizacáo 


t (5%) = 31 


0 T T T T T T T T T | Í 
0 2 4 6 8 10 12 14 16 18 20 


Figura 4.8 - Resposta ao degrau de um sistema sobreamortecido de segunda ordem. Os valores de fator de 
amortecimento estáo indicados nas curvas, considerando frequéncia natural unitária. 
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4.4.3 Sistemas criticamente amortecidos: [ = 1 


Para € = 1, o polinômio característico apresentado na Equação 4.2 possui um polo duplo em -o , 
e a resposta ao impulso se torna: 


h(t) = œ, t.e “ntu(t). (4.13) 
Tal resposta, mostrada nas Figuras 4.5 e 4.7, não apresenta oscilação, apresentando somente 
um período inicial de subida, seguido por uma descida assintótica que tende a zero. 


A resposta de um sistema criticamente amortecido ao degrau é dada por: 


s(D=[1-(1+0,t).en'Tu(t) (4.14) 


Tal resposta, ilustrada nas Figuras 4.6 e 4.8, corresponde a uma função crescente que tende 
assintoticamente ao valor unitário. 


Mais uma vez, o sobressinal da resposta ao degrau para sistemas criticamente amortecidos é 
nulo, pois a resposta nunca ultrapassa o valor máximo. Aproximações para os tempos de subida e de 
estabilização podem ser facilmente obtidas: 


» Tempo de subida 
t,(10-90%) = 3,47 
t, (595%) = 4,41 
» Tempo de estabilização 


t (5%) = 4,87 


4.5 Critério de estabilidade de Routh-Hurwitz 


O critério de estabilidade de Routh-Hurwitz constitui um método algébrico destinado a forne- 
cer informações sobre a estabilidade absoluta de um sistema linear invariante no tempo que possui 
um polinômio característico com coeficientes constantes. Tal critério avalia a possibilidade de uma 
das raízes desse polinômio se situar no semiplano lateral direito do plano complexo, ou seja, que 
possua parte real positiva, infringindo assim a condição primeira de estabilidade. 


O critério consiste na construção de uma tabela que estabelece uma condição necessária e 
suficiente de estabilidade, baseada em um arranjo triangular. Na tabela em questão, garante-se 
a estabilidade de um sistema se todos os elementos da primeira coluna do arranjo considerado 
forem positivos. Os passos do critério são descritos a seguir. 


a) Seja o polinômio característico: 


= n n-1 n=2 
H(s)=açs"+a  ,s”'+a, ,sP2+...+ajsta, 
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Amplie seus conhecimentos 


O critério de Routh-Hurwitz tem uma origem interessante. O algoritmo foi sugerido pelo matemático inglês Edward John 
Routh, Figura 4.9a, em 1876, com o nome de Teste de Routh, prestando-se a determinar se todas as raízes do polinô- 
mio característico de um sistema linear possuem parte real negativa. De modo independente, o matemático alemão Adolf 
Hurwitz, Figura 4.9b, propôs, em 1895, uma tabela de disposição dos coeficientes do polinômio característico em uma 
matriz quadrada chamada matriz de Hurwitz, mostrando que o polinômio é estável se e somente se a sequência de deter- 
minantes de suas submatrizes for positiva. 


2 c 
E £ 
E E 
E [o] 
8 (5) 
o o 
s = 
E E 
E Z 
z E 
5 E 

kej 
2 [=] 
2 e 
E 5 
o ho] 
2 3 
e 5 
2 < 


Figura 4.9a - Edward John Routh. Figura 4.9b - Adolf Hurwitz. 


Os dois algoritmos são absolutamente equivalentes, sendo o teste um modo elegante de se calcular os determinantes de 
Hurwitz. Um polinômio que satisfaça às condições de Routh-Hurwitz é, assim, chamado de polinômio de Hurwitz. 


b) Tomam-se os coeficientes de H(s), preenchendo-se a tabela que segue. 


a, dn2 n<4 


an- an-3 a, 5 
b; b, b; 
e; Cy C3 


c) Os elementos b; e c, são calculados de acordo com: 


bis An 1dn-2)7 An ân-2i-1 
i 
a 


n-l 


_ Dj An 217401 Pin 


C: 
1 b; 
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d) Completada a tabela, o número de mudanças de sinal na primeira coluna será o 
número de polos não negativos. 


O critério de Routh-Hurwitz é, assim, uma técnica eficiente no sentido de que permite a iden- 
tificação da estabilidade de um sistema sem o conhecimento exato da localização dos polos. 


Exercício resolvido 


2) Considerando-se o sistema com polinômio característico H(s) = sº +4s4 +25) +582 +3s+6, 
verifique sua estabilidade utilizando o critério de Routh-Hurwitz. 


Solução 


Monta-se a tabela de Routh conforme a seguir. 


[RS 0 
2 4 
E 
-3 6 

3 


Uma vez que há duas mudanças de sinal na primeira coluna, o sistema possui duas raízes no semi- 
plano lateral direito, no plano complexo, o que indica instabilidade. De fato, o sistema em malha 


fechada do exemplo possui raízes em s=-—3,7492, s = —0,6234+0,9054, 0,5205+1,0188, de 
modo que o último par de raízes, à direita do eixo imaginário, conduz o sistema à instabilidade. 


A maior vantagem do critério de Routh-Hurwitz é ser bastante funcional no projeto de con- 
troladores, já que a obtenção direta de valores que garantem a estabilidade não é trivial com o uso de 
cálculos simples. 


Exercício resolvido 


3) Considerando-se o sistema com ganho em malha aberta igual a: 


1000 


G()=————————— 
(s+2)s+3)(s+5) 


e realimentação com ganho unitário, verifique sua estabilidade utilizando o critério de 
Routh-Hurwitz. 
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Solução 
O sistema em malha fechada pode ser obtido como 


1000 

G(s) _ (s+2Ms+3X8+5) _ 1000 

1+G(s) 4, 1000 ~ $3 +1052 +31s+1030 ` 
(s+2)(s+3)(s+5) 


Humr (s) = 


A tabela de Routh é, assim: 


1 31 
10 1030 
2 0 
103 0 


== em) 


Uma vez que há duas mudangas de sinal na primeira coluna, o sistema possui duas raízes no 
semiplano lateral direito, no plano complexo, o que indica instabilidade. 


4.6 Análise dos resultados 


Agora que vocé já aprendeu um pouco mais sobre os sistemas de primeira e segunda ordem, 
algumas considerações podem ser levantadas. Em primeiro lugar, deve-se reforçar o fato de que, 
em sistemas de primeira ordem, quanto mais próximos os polos se situam do eixo imaginário, mais 
lentamente o sistema se estabiliza. Isto é reforçado no caso de sistemas subamortecidos, em que o 
tempo de estabilização é função da parte real dos polos do sistema. 


No que diz respeito a sistemas sobreamortecidos, a velocidade de estabilização depende uni- 
camente do polo mais próximo do eixo imaginário (dito polo dominante), sendo tanto mais lenta 
quanto mais próximo for esse polo do referido eixo. 


Outra conclusão relevante acerca dos sistemas apresentados é a de que sistemas com amorte- 
cimento crítico são os mais velozes em termos de convergência, posto que o polo mais próximo do 
eixo imaginário é o mais distante deste eixo, comparativamente aos outros tipos de sistemas. Isso é 
importante no projeto de sistemas mecânicos que necessitem de rápida estabilização, como sistemas 
que exijam rigidez mecânica (por exemplo, instrumentos de medição com pouca inércia). 


É possível também verificar que, em sistemas sobreamortecidos, polos mais próximos do eixo 
imaginário são considerados dominantes em termos de velocidade de resposta ao degrau unitário. 
Isso é também nítido quando se observam as respostas. 
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_ Vamos recapitular? 


Foram apresentados neste capítulo conceitos de estabilidade de sistemas lineares de primeira e 
segunda ordens, em termos de parâmetros considerados relevantes ao estudo das especificações requeri- 
das para sistemas de controle. 


O capítulo abordou, entre outros elementos, os conceitos de tempo de subida, tempo de pico, 
tempo de estabilização e sobressinal, que constituem elementos caracterizantes de sistemas lineares 
quando sujeitos a entradas impulsivas e do tipo degrau. Analisou-se também um critério de verificação 
de estabilidade de sistemas, o critério de Routh-Hurwitz, que dispensa para tal o conhecimento exato da 
localização dos polos do sistema em análise. 


O Agora é com você! 


1) Avalie a estabilidade dos sistemas representados pelas funções de transferência abaixo: 


3 
H(s) === 
6) s? +5s+16 
4 
H(s) = ==—— 
(5) s? +4s+10 
3 
H(s) = 
(9) s2 +1,8s +2,2 
: 2) Determine para o sistema descrito por E os valores de: 
: 2s? +10s+14 


a) tempo de pico; b) sobressinal percentual; c) tempo de acomodação. 


: 3) Sob que condições um sistema de segunda ordem tem sua saída aproximada à de um 
sistema de primeira ordem? Justifique sua resposta. 


K 
8 
+ SA 


; 4) Para o sistema em malha aberta G(s) = > 
: (s+2)(s+3) 


referente a um sistema controlado em modalidade proporcional integral (a ser visto 
nos próximos capítulos), determine os valores das constantes envolvidas, utilizando 
o critério de estabilidade de Routh-Hurwitz, para que o sistema em malha fechada 
seja estável. 
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Técnica do Lugar 
das Raízes (LGR) 


Para começar 


Este capítulo tem por objetivo apresentar a técnica do Lugar das Raízes (LGR), que permite a aná- 
lise de comportamento de um sistema em malha fechada, visando a observar todos os seus possíveis 
valores possíveis de ganho a partir da movimentação de seus polos. 


As informações contidas neste capítulo constituem uma relevante base para o projeto de contro- 
ladores a partir da localização de raízes do sistema, posto que a aferição de ganho de acordo com tais 
características é de fácil acesso ao projetista. 


5.1 Introdução ao diagrama LGR 


O diagrama do Lugar das Raízes (LGR, também chamado de diagrama de Evans ou root locus) 
consiste em um conjunto de curvas no plano complexo s, em que tais curvas representam as posi- 
ções admissíveis para os polos de malha fechada de um dado sistema, na condição de que seu ganho 


varie de zero a infinito. 


Considere-se a Figura 5.1, indicativa de um sistema realimentado. 
O ganho em malha fechada de tal sistema pode ser escrito como: 


C(s)_ G(s) (5.1) 
R(s) 1+G(S)H(s) ` 


Gyme(s) = 


R(s) C(s) 
“O o 


H(s) 


Figura 5.1 - Diagrama de blocos de um sistema em malha fechada. 


Os polos de malha fechada são as raízes do polinômio característico 1+G(s)H(s), ou seja, 


para sua determinação é necessário que G(s)H(s) =—1, o que leva às condições: 


|G(s)H(s)=1 (5.2) 


G(s)H(s) =+180(2k+1), k=0,1,2... (5.3) 


Deste modo, o objetivo da técnica é determinar os polos de malha fechada em questáo, uma vez 
que estes caracterizam a resposta do sistema analisado. Neste caso, o parámetro K se torna explícito 
no sistema em malha aberta, como será visto nos próximos exemplos, sendo relevante ao projetista. 


5.2 Projeto de controladores por lugar das raízes - 
passos de elaboracáo 


Para determinação do Lugar das Raízes, o processo não é trivial, pois envolve o conhecimento 
apurado de matemática aplicada ao processo. Para levantamento das curvas do Lugar das Raízes 
associado a um sistema em particular, devem ser executados passos específicos, descritos a seguir. 


1. Escrever o polinômio característico de modo que o pará- 1 
+ H(s) =1+K.P 
metro de interesse (K) apareça claramente. Exemplo: NH) (9) 


Tietz 


2. Fatorar o polinômio P(s) em termos dos np polos e 1+G(s)H(s)=1+K El — 


P 
ny zeros. 
Z I I (s+ pj) 
FI 


3. Assinalar os polos e zeros de malha aberta no plano s 
com os símbolos correspondentes. 


4. Assinalar os segmentos do eixo real que são LGR. O LGR se situa à esquerda de um número ímpar de polos e zeros 


X — polos, O — zeros (o LGR começa nos polos e termina nos zeros) 
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LS=np, quando np >n7 


5. Determinar o número de lugares separados, LS (seg- Np = número de polos finitos 
mentos de curva que compóem o LGR). 
nz = número de zeros finitos 


E para É eisimétrico jcomirelação aofeixo reali(eixo hori- Basta desenhar a parte acima do eixo real e depois espelhar o esboço. 


NO) IE) 
A" 


7. (np —nz) segmentos de um LGR prosseguem em dire- omiy 
ção aos zeros infinitos, ao longo de assíntotas centrali- 
zadas em o, e com ângulos ba - da 


_ q+!) 
gamy 


“180º, =0,1,2,....(np —nz —1) 


a) Fazer K =P(s); 
8. Determinar o ponto de saída (se existir) sobre o eixo aña ©) 


real. b) Determinar as raízes de dP(s) _ 0 
ds 
9. Utilizando o critério de Routh-Hurwitz, determinar 


o ponto no qual o eixo imaginário é cruzado (se isso | Ver critério de estabilidade de Routh-Hurwitz 
ocorrer). 


Ángulo de partida = 180° D'0,+)'9; 


Ángulo de chegada = 180°- q, +70; 
10. Usando a condição de ângulo, determinar o ángulo de 
partida para os polos complexos e o ángulo de chegada | em que: 
para os zeros complexos. 
0; = ángulos de vetores partindo dos demais polos até o polo em questão 


Q; = ângulos de vetores partindo dos demais zeros até o polo em questão 


Exercício Resolvido 


1) Esboce o LGR para o sistema em malha aberta G(s) = RE 
s(s+4) 


Solução 


1. Escrever o polinômio característico de modo que o parâme- 14+G(s)H(s)=1+ K(s+2) 


tro de interesse (K) apareça claramente s2+4s 


KUI 0) > s+2 
s? +4s s(s+4) 


2. Fatorar o polinômio P(s) em termos dos np polos e nz | 1+G(s)H(s)=1+ 
zeros. 
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3. Assinalar os polos e zeros de malha aberta no plano s com 
os símbolos correspondentes. 


Parte imaginária 


| L I 1 I 
35 3 -25 2 15 


Parte real 


Lugar Geométrico das Raízes 


2 polos e zeros 
(nº par) 


-4 i -3 2 


3 polos e zeros 
(nº impar) 


4. Assinalar os segmentos do eixo real que sáo LGR. 


1 polos e zeros 
(nº ímpar) 


LGR final 


Root Locus 


ginário 


o 


Eixo ima 


14 12 10 8 6 4 —2 0 2 
Eixo real 
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(9) Exercício Resolvido 


2)  Esboce o LGR para o sistema em malha aberta H(s) = e o , com ganho de 
(s+2)(s+ 4) 


realimentação igual a A 
bs s(s+4) ` 
Solução 


1. Escrever o polinômio característico de modo que o parâme- pá K(s+1) 
tro de interesse (K) apareça claramente EG do s4 +1083 +3282 +32s 


2. Fatorar o polinômio P(s) em termos dos np polos e nz 2 s+1 
s(s+2)s +4)? 


P(s) 


zeros. 


Polo com 
multiplicidade 2 


3. Assinalar os polos e zeros de malha aberta no plano s com 
os símbolos correspondentes. 


Trecho entre 
2 polos 
Polo com Total de 
multiplicidade 2 2 polos e zeros 


(nº par) 


4. Assinalar os segmentos do eixo real que são LGR. 
-5 4 -3 —2 —1 0 


Total de 
3 polos e zeros 
te ímpar) Total de 
1 polos e zeros 
(n° impar) 


sl di AO 


6. O LGR é simétrico com relacáo ao eixo real 
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7. (np —nz) segmentos de um LGR prosseguem em direção 
aos zeros infinitos, ao longo de assíntotas centralizadas em 
Ga e com ángulos ba - 


o, -È 
E Mp Mz 
ESA 
4-1 3 
q+) 
Epi 


180º, q=0,1,2,...,(np -nz —1) 


60º, q=0 
180º, q=1 
300º, q=2 


dP(s) _ 2 
8. Determinar o ponto de saída (se existir) sobre o eixo real. ds 0 > s=-2,5994 


Fazer K = P(s); (Ponto de saída sobre Re) 

Determinar as raízes de dP(s) _ 
ds 

K(s+1) 
s4 +1083 + 328? + 325 
s4 +1083 + 328? + 32s 

a s+1 

dP(s) _ 354 + 2483 +6252 + 64s +32 E 

FARS (s+1? E 

s=-2,5994 


0 
1+K.P(s) =1+ 
P(s)=K = 


0 


LGR final 


0.8 i l ioa: lecie 


Eixo imaginário 


Eixo real 
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Fique de olho! 


O LGR é sempre simétrico em relação ao eixo real, já que os polos de um sistema linear, quando não são reais, sempre 
são complexos conjugados. 


5.3 Lugar das raízes para funções de 
transferência típicas 


O Quadro 5.1 a seguir sumariza diferentes lugares das raízes para algumas das funções de 
transferência usuais no plano complexo. 


Quadro 5.1 - Lugar das raízes para funções de transferência usuais no plano complexo. 


KG(s)H(s) 


E Sos 
(ts+D(t,s+1) 


RO A 
(us +1D(r)s +1)(135+1) 
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E c 
s 

-5 

s(qs+1) 


do. 
s(1,s+1)(t,s +1) 
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K(st, +1) 
s(ņqs +1)(t,s +1) 


Polo com 
multiplicidade 2 


Polo com 
multiplicidade 2 


-E 
s? (qıs +1) 
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K(st, +1) 
n G, 
s? (qs +1) 


>7 


Polo com 
multiplicidade 3 


Polo com 
multiplicidade 3 


K(st, +1) 
s3 
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Polo com 
multiplicidade 3 


K(st, +1)(61, +1) 
3 
s 


Polo com 
multiplicidade 2 


K(st, +1) 


s? (st; +1)(st, +1) 


K(st, +1)(s7, +1) 
(st; + 1)(s7, + 1)(s7 +1)(st4 +1) 
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Fique de olho! 


Os pontos de partida do LGR (pontos em que K=0) coincidem com os polos da função de transferência em 


malha aberta. 
( Vamos recapitular? 


Foram apresentados neste capítulo conceitos para a construção do lugar geométrico das raí- 
zes, técnica útil na observação do comportamento de um sistema a entradas conhecidas. Por meio da 
variação de ganho em malha aberta, é possível observar os valores de resposta em malha fechada, o 
que serve de suporte ao projeto de controladores. Igualmente, o LGR oferece auxílio à avaliação das 
características da resposta no tempo do sistema em malha fechada, como função da variação de parâ- 
metros conhecidos. 


O capítulo reforçou, entre outros elementos, a relevância do conhecimento da localização das raí- 
zes do polinômio característico e suas condições de obtenção a partir do manuseio do LGR. 


: Agora é com você! 


1) Esboce o lugar das raízes para os sistemas em malha aberta apresentados a seguir. 


a) G(s) = += 

b) G(s) = Ds 

c) G(s)= E E 
d) G(s)= EA 


É possível, para tais sistemas, designar diretamente no gráfico a posigáo em que as cur- 
vas “explodem”? Justifique sua resposta. 


: 2) Para o lugar das raízes ilustrado a seguir, o que se pode afirmar a respeito do sistema 
em malha fechada considerado? 
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K + da) 
a PA 


3) Para o sistema em malha aberta G(s) = 
(s+2)(s+3) 


, apresentado no capítulo anterior: 


» atribua valores para as constantes que assegurem estabilidade, e esboce o lugar 
das raízes do sistema arbitrado. 


4) O controle da pressão de sucção de um compressor é apresentado na figura a seguir. 


setpoint 


Considerando-se constante a vazão de entrada no reservatório (f; ), pode-se deduzir o 
seguinte modelo dinámico que relaciona a vazão no compressor f, (em unidades físicas 
adequadas) à pressão de sucção ( Pç, expressa em psi): 


—0.5 


E: + 
s(8) 7,5s+1 


F (s) 


Além disso, a resposta da vazão no compressor em relação ao sinal de comando M(s) 
(em %) é dada por: 


0,36 
E 


O transmissor de pressáo possui ganho unitário e constante de tempo de 1,2s. 
Com base nas informações dadas, responda aos itens a seguir. 

a) O controlador deverá possuir ação direta ou realimentada? Explique. 

b) Faça um diagrama de blocos do sistema em malha fechada. 


c) Obtenha a função de transferência em malha fechada. 
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Ações de Controle 


Para começar 


Este capítulo tem por objetivo apresentar as ações de controle mais conhecidas, de acordo com as 
especificações de desempenho desejadas. São abordados conceitos de controle e compensação, à luz das 
modalidades de controle proporcional, controle derivativo e controle integral, além das técnicas de com- 
pensação baseadas em avanço e atraso de fase. 


As informações contidas neste capítulo constituem o cerne da tarefa do projetista de controle de sis- 
temas lineares, sintetizando toda a teoria apresentada ao longo desta obra. No caso, o uso de uma estratégia 
sobre outra pode revelar vantagens e desvantagens que não devem escapar ao entendimento do projetista. 


6.1 Estratégias de controle 


Como expresso no Capítulo 1, um controlador tem por finalidade manipular a relação entre 
sinal de saída e sinal de entrada de um sistema, atuando sobre seus parâmetros a fim de se atingirem 
especificações de desempenho. Por meio da manipulação de algumas variáveis, controlam-se outras, 
no que o controlador pode ser representado sob diversos enfoques, advindo daí a necessidade de o 
projetista saber qual modalidade de controle utilizar. Nesse sentido, um controlador pode atuar em 
malha aberta ou malha fechada (realimentado). 


Dentro da classe de controladores realimentados, destacam-se duas abordagens: as que ope- 
ram sobre o ganho (módulo) dos sinais envolvidos e as que operam sobre a fase destes. Faz-se aqui 
uma pormenorização dessas abordagens, para estudo. 


6.2 Controle on-off 


Nesta estratégia de controle, o elemento de atuação possui apenas dois estados. Sua aplicação 
se justifica pela facilidade de construção e baixo custo. Dependendo do sinal do erro, a entrada é 
estabelecida na posição máxima ou mínima. 


O segundo termo é mínimo e igual a zero ou o oposto do máximo. Tais controladores são 
geralmente implementados com o uso de válvulas solenoides. 


6.3 Controle proporcional (P) 


A ação de controle proporcional estabelece uma relação direta entre saída e entrada, dada pelo 
ganho proporcional K , mostrado na Figura 6.1. O erro de estado estacionário é denominado egg, 
e é inversamente proporcional ao ganho. Todavia, a diminuição do erro conduz a um maior tempo 
de estabilização, podendo até conduzir à instabilidade. Por outro lado, a redução do ganho torna o 
sistema menos afeito à instabilidade, reduzindo-se também o tempo de estabilização, mas elevando 
o erro de regime. Um sistema de controle do tipo proporcional não acrescenta zeros nem polos ao 
sistema em estudo, constituindo apenas um ajuste de ganho. 


u(t) = Ke(t) = K[r(t) — y(t)] 


(6.1) 
U(s) = KE(s) 
C(s) 
Figura 6.1 - Diagrama de blocos de um controlador proporcional atuando em malha fechada. 
Para uma entrada do tipo degrau, o erro em estado estacionário se resume a: 
Ess =7 > (6.2) 
SS 1+K 


o que nem sempre é viável. 


6.4 Controle proporcional + derivativo (PD) 


A ação derivativa, associada à ação proporcional, provoca o acréscimo de um zero ao sistema, 
o que é benéfico na resposta transitória, além de fornecer maior estabilidade ao sistema e reduzir o 
tempo de estabilização, já que se insere um zero em 
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Por outro lado, eleva-se o tempo de subida, além de não corrigir o erro em regime permanente. 


sn) 
dt 


u(t) = Kp (eo +Ta 


U(s) = (Kp +K ¿s)E(s) (6.3) 
Ki =Kpta 


A ação derivativa, de fato, introduz um efeito de antecipação na resposta do sistema, fazendo 
com que ele responda não somente à magnitude do sinal de erro, mas também à sua tendência para 
tempos posteriores, agindo preventivamente, Figura 6.2. Uma vez que a derivada de um sinal repre- 
senta sua inclinação da função, o controle derivativo é essencialmente um tipo que antecipa a ten- 
dência desse sinal, predizendo um sobre ou subsinal à frente no tempo, perfazendo as apropriadas 
correções antes daquele sobre ou subsinal acontecer. 


R(s) C(s) 
+ Ko + Kys SA 


Figura 6.2 - Diagrama de blocos de um controlador proporcional + derivativo em malha fechada. 


6.5 Controle proporcional + integral (PI) 


A ação de controle proporcional associada à integral corresponde ao estabelecimento de uma 
taxa de variação do sinal de saída com relação ao sinal de entrada na forma: 


t 
u()=Kp| e) ++ fedir] > u()- CEO > E= (6.4) 
S 


o 1 


Com a ação integral, atua-se de modo efetivo na resposta em regime permanente, o que dimi- 
nui o erro de estado estacionário. Todavia, traz-se prejuízo ao regime transitório, pois são acrescidas 
raízes ao sistema que tendem a desestabilizá-lo (um polo em p=0 e um zero em z=-—+ 

K ), o que 


também acarreta aumento do tempo de estabilizacáo. 


Ao se condensar a ação proporcional à integral, é possível um melhor desempenho na resposta 
transitória em razão da ação proporcional, ao passo que a ação integral age sobre o erro em regime 
permanente, em que a ação proporcional falha, como mostra a Figura 6.3. 


R(s) C(s) 


K 
AL 
i Kg s ts +1 


Figura 6.3 - Diagrama de blocos de um controlador proporcional + integral atuando em malha fechada. 


Ações de Controle 


Para uma entrada do tipo degrau, o erro em estado estacionário se resume a: 


1 
ess E =0, (6.5) 


Tal resultado é esperado, pois a finalidade da ação de controle integral visa justamente a elimi- 
nar o erro em regime permanente. 


6.6 Controle proporcional + integral + derivativo (PID) 


A ação proporcional aliada às ações integral e derivativa permite efeitos tanto no regime transitório 
quanto no regime permanente, prestando-se assim a solucionar problemas nos dois regimes, Figura 6.4. 


t 
u(t) =K | e(t) de f e(T)dT+ T4 
To 


(Kps+K, +K s?)E(s) 


de(t) 


U(s) = (6.6) 


K 
K= 2, Ka = Kpta 
Ti 


K 
=L E. 
+ Kp + Kgs + A ma! a 


Figura 6.4 - Diagrama de blocos de um controlador proporcional + integral + derivativo atuando em malha fechada. 


Boa parte dos controladores industriais atuais faz uso da estratégia PID. Isso se deve ao fato de 
que o desempenho de tais controladores é robusto sobre uma ampla faixa de condições de operação, 
além de sua simplicidade de implementação. 


A Tabela 6.1 sumariza as características principais das ações proporcional, integral e deriva- 
tiva, atuantes simultaneamente em um controlador PID. 


Tabela 6.1 - Especificações de desempenho para cada modalidade de controle realimentado 


Ganho Tempo de subida Sobressinal Tempo de estabilização Erro de regime permanente 


so a C 


gi EEE COC ii TEN 


O efeito resultante de fato da operação de um controlador PID depende, entretanto, da preci- 
são do projeto, podendo superar os efeitos citados na Tabela 6.1 como referência. 
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6.7 Controle por avanço de fase (LEAD) 


A ação de controle por avanço de fase tem por fim compensar um eventual atraso de fase 
surgido dentro do sistema a ser controlado, por meio do redesenho do lugar das raízes a fim de se 
obterem os polos adequados em malha fechada. O uso da compensação LEAD provoca aumento 
no amortecimento do sistema, reduzindo-se os intervalos de subida e de estabilização, sem alteração no 
erro de regime permanente. 


A função de transferência correspondente é: 


ts+1 
Ge(s)=Kc 


, O<a<l (6.7) 
ots+1 


Em suma, a ação de controle LEAD introduz um zero e um polo no sistema, melhorando o 


desempenho em regime transitório (tal como no controlador PD), à custa de um avanço na fase total 
do sistema. 


6.8 Controle por atraso de fase (LAG) 


A ação de controle por atraso de fase (LAG) traz melhora quanto ao erro em regime estacioná- 
rio, acarretando, todavia, uma resposta mais lenta - maiores intervalos de subida e de estabilização. 
Por outro lado, o sistema se torna mais robusto à ação de ruídos espúrios. 


A função de transferência correspondente é: 


ts+1 
Go(s) = Ko Pesa? B >1 (6.8) 


Em suma, a ação de controle LAG introduz um zero e um polo no sistema, melhorando o 


desempenho em regime permanente (tal como no controlador PI), à custa de um atraso na fase total 
do sistema. 


6.9 Controle por avanço-atraso de fase (LEAD/LAG) 


A ação de controle por avanço-atraso de fase (LEAD/LAG) tem por finalidade a obtenção da 
resposta rápida oferecida por compensadores LEAD, associada a uma diminuição do erro em regime 
permanente, garantida por compensadores LAG. 


A função de transferência correspondente é: 


+1 +1 
PE japa pod, Las (6.9) 
arts +1ft,5+1 v H 


Ge(s) =Kc 


Em suma, a ação de controle LEAD/LAG introduz dois zeros e dois polos no sistema, melho- 
rando o desempenho tanto em regime transitório como em regime permanente (tal como no con- 
trolador PID). 


Ações de Controle 


_ Vamos recapitular? 


Foram apresentados neste capítulo conceitos referentes às modalidades de controladores estuda- 
dos em sistemas de controle: controlador P, controlador PI, controlador PD, controlador PID, compen- 
sador LEAD, compensador LAG, compensador LEAD/LAG. 


O capítulo abordou, entre outros elementos, as características de cada tipo de controlador, listan- 
do as vantagens e desvantagens de cada um. Tais considerações são relevantes no sentido de se arbitrar, a 
partir das especificações de desempenho desejadas, qual estratégia de controle utilizar. 


O Agora é com você! 


1) Dadas as funções de transferência dos seguintes sistemas (malha direta e malha rever- 
sa), traçar o LGR do sistema sem o controlador, testando se o ponto dado pertence ao 
LGR do sistema. A seguir, traçar o LGR com o sistema em série com um controlador 
PD com ganhos Kp = 0,4 eK, = 0,1, testando se o ponto dado pertence ao LGR. 


1 1 
G(s) = —— H(s) = 5 =-—(0,3044+3,438i 
(s) úl (s) z s i +3,438i 
: 2) Considere o sistema em malha aberta G(s) = ml i 
: (s+2)(s+5) 


Verifique se o ponto s =-2+ j3 pertence ao LGR do sistema em questão para algum 
valor de K . Considere realimentação com ganho unitário. 


: 3) Para os sistemas mostrados abaixo, aponte, qualitativamente, as diferenças no tra- 


çado do LGR de cada um. 
R(s) a c(s) 
Ei s(s+8) á 
R(s) = cl) 
+ s(s+8) 
1 
s(s + 4) 
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Projeto de 
Controladores 
por Meio do 
Lugar das Raízes 


Para começar 


Este capítulo tem por objetivo apresentar os conceitos de projeto de controladores com base no 
traçado do Lugar das Raízes. O capítulo traz o detalhamento de projeto dos controladores já descritos 
no capítulo anterior, com uma nova forma de sua visualização por meio da aplicação da técnica do Lugar 
das Raízes. 


As informações contidas neste capítulo trazem um complemento aos estudos apresentados no 
capítulo anterior, de forma a aperfeiçoar as regras de projeto já enunciadas. 


7.1 Parâmetros do controlador 


No projeto de controladores, são usualmente empregados parâmetros referentes a sistemas de 
segunda ordem, como já explicitado nos capítulos anteriores. Nesse sentido, são relevantes para pro- 
jeto as especificações de fator de amortecimento £ e frequência natural ©, . 


O polinômio característico de um sistema de segunda ordem é dado por 


P(s) =s? +2f0,5+0,?, (7.10) 
com polos em 


s=—Lo, +jo, y1-8? q, (7.2) 


O projeto de controladores usualmente se restringe a um método em força bruta, em que se 
buscam atingir as especificações de desempenho requeridas em malha fechada iterativamente até 
que sejam obtidas, por meio da modificação dos parâmetros livres para ajuste ou mesmo da configu- 
ração do sistema. 


Quando se deseja alterar o desempenho do sistema em regime transitório, o controlador apli- 
cado deve inserir raízes no sistema global de modo que o lugar das raízes cruze o ponto especificado 
de desempenho. Em regime permanente, o controlador projetado deve operar sobre o ganho do sis- 
tema, promovendo poucas modificações no lugar das raízes original. 


7.2 Projeto de controladores PD com base no LGR 


Recordando o capítulo anterior, a função de transferência de um controlador PD é dada por: 


GeO- ES = Ke +Kys=K¿(s+2) (7.3) 


Apontam-se os seguintes passos de projeto para controladores PD a partir da observação do LGR: 


U(s 
s 


a) A partir das especificações de desempenho requeridas, obter a localização desejada dos 
polos dominantes de malha fechada. 


b) Testar a viabilidade de uso de um controlador P puro. 


c) Caso seja preciso utilizar um controlador PD, posicionar o novo zero de modo a satisfazer 
a condição de ângulo do LGR. 


d) Obter o ganho total requerido a partir da condição de módulo. 
e) Determinar o valor da constante de erro estacionário. 
f) Caso tal constante não seja viável, alterar o tipo de controlador. 


g) Simular o sistema com o controlador e observar o comportamento da resposta, ajustando 
de modo preciso os parâmetros Kọ e z. 


Em termos práticos, pode-se utilizar a regra empírica de projeto: 


_U(s) + K ¿s 
IN y 3<a<l0 (7.4) 
1+= 
a 
7.3 Projeto de controladores Pl com base no LGR 
A funcáo de transferéncia de um controlador PI é dada por: 
U(s) K. S+Z 
Go(s)= =K,+2=K Te, 
c(s) mo a ta (7.5) 


Apontam-se os seguintes passos de projeto para controladores PI a partir da observação do LGR: 


a) Localizar o polo situado na origem. 
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b) Determinar a localização do zero, de maneira a satisfazer a condição de ângulo do LGR. 
c) Determinar o ganho total desejado, por meio da aplicação da condição de módulo. 
d) Simular a resposta do sistema em malha fechada com o controlador atuando. 


e) Caso o desempenho do controlador não seja satisfatório, realizar um ajuste fino dos parâ- 
metros do controlador (Kọ e z). 


7.4 Projeto de controladores PID com base no LGR 


A função de transferência de um controlador PID é dada por: 


E: 
GeO py Kpt +K = 


s?K4a+sKp +K; E (s+z¡Xs+Z)) 


; 3 A (7.6) 


Apontam-se os seguintes passos de projeto para controladores PID a partir da observação do LGR: 


a) A partir das especificações de desempenho, obter a localização desejada dos polos dominan- 
tes de malha fechada. 


b) Observar se é possível a obtenção de tais especificações com o uso de controladores P, PI 
ou PD. 


c) Caso seja necessária de fato a utilização do PID, determinar a localização do polo na ori- 
gem e dos zeros tal que a condição de ângulo do LGR seja satisfeita. 


d) Determinar o ganho total do sistema, a partir da condição de módulo do LGR. 


e) Simular a resposta com o uso do controlador projetado. Caso o desempenho do controla- 
dor não seja satisfatório, realizar um ajuste fino dos parâmetros do controlador (Kç, z; 
ez). 
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7.5 Projeto de controladores por avanço 
de fase com base no LGR 


A função de transferência de um controlador por avanço de fase é dada por: 


1 
s+— 
atj- —L E penal (7.7) 
E(s) PP s+p 
aT 


Apontam-se os seguintes passos de projeto para controladores por avanço de fase a partir da 
observação do LGR: 


a) A partir das especificações de desempenho, obter a localização desejada dos polos domi- 
nantes de malha fechada. 


b) Observar se é possível a obtenção de tais especificações com o uso de um controlador do 
tipo P (proporcional). 
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c) Caso seja necessária de fato a utilização do avanço de fase, determinar a localização do zero 
do controlador, bem como posicionar o polo do controlador de tal forma que a condição 
de ângulo do LGR seja satisfeita. 


d) Determinar o ganho total do sistema, a partir da condição de módulo do LGR. 
e) Determinar o valor da constante de erro estacionário. 
f) Sea constante não for adequada, realizar nova escolha de controlador e retornar ao passo (c). 


g) Simular a resposta com o uso do controlador projetado. Caso o desempenho do controla- 
dor não seja satisfatório, realizar um ajuste fino dos parámetros do controlador (Kọ e z; ). 


7.6 Projeto de controladores por atraso 
de fase com base no LGR 


A função de transferência de um controlador por atraso de fase é dada por: 


1 
= == poa TA B>1 (7.8) 


Apontam-se os seguintes passos de projeto para controladores por atraso de fase a partir da 
observação do LGR: 


a) Obter o ganho de malha aberta por meio do uso da condição de módulo. 
b) Determinar o valor da constante de erro estacionário. 


c) Determinar o quanto deve ser elevada tal constante para que se atinjam as especificações 
de desempenho requeridas. 


d) Obter o polo e o zero do controlador que produzem o devido aumento, sem alteração sig- 
nificativa do LGR original. 


e) Simular a resposta com o uso do controlador projetado. 


f) Se preciso for, remodelar as raízes do passo (d) até que se atinjam as especificações de 
desempenho requeridas. 


7.7 Projeto de controladores por avanço-atraso 
de fase com base no LGR 


A função de transferência de um controlador por avanço-atraso de fase é dada por: 


R so sto 
Go(s) -YD cj bj b (7.9) 
c c 
E(s) A 
T BT, 
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Apontam-se os seguintes passos de projeto para controladores por avanço-atraso de fase a par- 
tir da observação do LGR: 


a) A partir das especificações de desempenho, obter a localização desejada dos polos domi- 
nantes de malha fechada. 


b) Determinar a contribuição de ângulo necessária para a rede LEAD. 

c) Determinar as constantes T, e y por meio da condição de ángulo do LGR. 
d) Determinar a constante K¿ por meio da condição de módulo do LGR. 

e) Determinar a constante B a partir das especificações de erro. 


f) Determinar a constante T, de modo que a contribuição de ángulo da seção LAG se situe 
entre -5 e 0º (regra prática). 


g) Simular a resposta com o uso do controlador projetado, procedendo ao ajuste fino dos 
parâmetros do controlador, se preciso for. 


( Vamos recapitular? 


Foram apresentados neste capítulo conceitos referentes às modalidades de controladores estuda- 
dos em sistemas de controle: controlador P, controlador PI, controlador PD, controlador PID, compensa- 
dor LEAD, compensador LAG, compensador LEAD/LAG, desta vez sob o ponto de vista da construção 
de traçados de lugar das raízes. 


O capítulo abordou, entre outros elementos, técnicas alternativas ao projeto de tais controladores, 
enunciando um conjunto de passos para tal, de modo a desfrutar da funcionalidade na construção do 
LGR. 
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O Agora é com você! 


: 1) Para os sistemas dados a seguir, projetar controladores que atendam às especificações 
: de desempenho requeridas. 
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A Es 2000(s? —1,14) 


Projetar um controlador PD para que SS=0,8 e q, =0,4rad/s em malha 
fechada. 


o 100 
—s(s+8Xs+24) 


Projetar um controlador em atraso para que o fator de amortecimento dos 


b) G(s) 


polos dominantes se situe em € = 0,707. 


: 2) Para o sistema G(s)= SS i 
: s2 +6s+9 


obtenha por simulação a resposta do sistema após a atuação de um controlador PI, com 
os parâmetros: 


Kc=1, z=0,01 
Ko=1 z=0,3 
Kc =3, z=0,3 
Kc =9, z=0,9 
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Projeto de 
Controladores 
no Domínio da 
Frequência - 
O Diagrama 
de Bode 


Para começar 


Este capítulo tem por objetivo apresentar os conceitos de projeto de controladores no domínio 
da frequência, que avalia a resposta de sistemas controlados no que se refere ao seu comportamento de 
módulo e fase a partir de respostas obtidas a entradas conhecidas, com ênfase no uso de Diagramas 
de Bode. 


As informações contidas neste capítulo fornecem um subsídio diferenciado ao entendimento de 
sistemas de controle, apresentando uma abordagem (domínio) diferente na obtenção de especificações 
desejadas de projeto. 


8.1 Introdução ao projeto de controladores 


Quando se faz uso do termo “resposta em frequência” em Sistemas de Controle, busca-se 
avaliar a resposta em regime permanente de um sistema a uma entrada tipicamente senoidal, em 
que o sinal de entrada tem sua frequência natural variada a fim de se observar as características de 
resposta daquele sistema. Nesse sentido, as conclusões retiradas de tal análise são um tanto distin- 
tas das obtidas quando se faz uso de técnicas como a do Lugar das Raízes. A maior vantagem de 
métodos baseados na análise em frequência reside no fato de que os dados obtidos da simulação 
de plantas sob tal enfoque dispensam o uso de modelos matemáticos levantados para caracterizar 
tais processos, conquanto sejam úteis apenas na análise temporal desses processos. 


Há uma variedade de métodos de controle de processos baseados na análise da resposta 
em frequência, concentrados na primeira metade do século passado, como os desenvolvidos por 
Nyquist, Bode e outros estudiosos da Teoria de Controle. 


Amplie seus conhecimentos 


Hendrik Wade Bode 


O engenheiro Hendrik Wade Bode (1905-1982), funcioná- 
rio da Bell Labs na década de 1930, trouxe valiosas con- 
tribuições à Teoria de Controle. De suas ideias, uma das 
mais relevantes para a análise da resposta de sistemas no 
domínio da frequência é o chamado Diagrama de Bode, con- 
junto de gráficos de módulo e de fase de um sistema sujeito 
a entradas com frequências variadas dentro de uma faixa de 
trabalho específica. 


Wikipedia/Commons 


Figura 8.1 - O engenheiro Hendrik Bode, idealizador dos 
Diagramas de Bode. 


Como vantagem maior das técnicas de controle baseadas no domínio da frequência, aponta-se 
a possibilidade de que funções de transferência complexas possam ser experimentalmente modela- 
das, pelo uso de simulações. Não obstante, torna-se possível um melhor tratamento de condições 
em que se faça presente a contaminação por ruído, posto que a caracterização em frequência admite 
tal possibilidade. 


A próxima seção destaca a implementação dos Diagramas de Bode como estratégia mais usual 
de controle baseada no domínio da frequência. 


8.2 Diagramas de Bode 


O Diagrama de Bode consiste em uma construção gráfica com duas curvas distintas: a pri- 
meira representa o logaritmo, na base 10, do módulo de uma função de transferência representada 
no domínio da frequência, ao passo que a segunda representa o ângulo de fase dessa mesma função 
de transferência. O gráfico de módulo, no Diagrama de Bode, é expresso em dB (decibel), posto que 
se traça uma curva logarítmica 20log|G(jw)|, em que G(s)é a função de transferência analisada, no 
domínio da transformada de Laplace. 


Fique de olho! 


A representação gráfica dos Diagramas de Bode, tanto para a curva de módulo quanto para a curva de fase, requer a 
representação do eixo das frequências em escala logarítmica. Assim, se você resolver esboçar os diagramas em papel, 
utilize papel monolog para os gráficos de fase, e papel dilog (ou log-log) para os gráficos de módulo. 
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A utilização de Diagramas de Bode encontra como principal vantagem o fato de que eventuais 
produtos aritméticos entre módulos podem ser convertidos em somas, em razão das propriedades 
dos logaritmos. Além disso, é possível representar as curvas por meio de aproximações assintóticas, 
bastante úteis à compreensão do comportamento dos sistemas analisados conforme a frequência. 


A construção de Diagramas de Bode segue como premissa a separação da função de transferên- 
cia analisada em fatores, que permitem a construção dos gráficos de forma simples e gradual. A com- 
binação desses fatores por soma simples permite a representação de funções mais complexas, já que 
somar tais fatores corresponde, em termos logarítmicos, ao produto dos módulos. Vejamos, nas pró- 
ximas subseções, quais os fatores que podem ser representados na construção de Diagramas de Bode. 


8.2.1 Fator de ganho K 


A representação de valores constantes na escala logarítmica consiste em representar valores 
maiores que a unidade como positivos, e valores menores que a unidade como negativos. Assim, 
uma constante K é representada na curva de módulo do Diagrama de Bode como uma reta horizon- 
tal com 20logKdB. Valores constantes possuem valores de fase nulos. Desse modo, uma constante 
apenas eleva ou diminui o nível da resposta de módulo de um valor de 20logKdB, não promovendo 
nenhuma alteração na fase. A Figura 8.2 representa a curva de módulo para uma constante K, con- 
forme a construção definida pelo Diagrama de Bode. 


Módulo (dB) 


0,001 0,01 0,1 di 10 


Fase (graus) 


Li Lu aaa) 
0,001 0,01 


Oo 10 
Log (frequência) 


Figura 8.2 - Diagramas de Bode para um fator de ganho K = 10. 
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8.2.2 Fatores derivativo e integral (jw)*? 


O fator 20logljw| é igual a 20log(w)dB, o que corresponde a uma reta inclinada em 20 dB/ 
década na curva de módulo do Diagrama de Bode, bem como a um valor constante positivo na 
curva de fase do diagrama. Do mesmo modo, o fator 20log|(jw) 'lequivale ao termo -20log(o), 
o que implica raciocinar, para a curva de módulo do Diagrama de Bode, uma reta inclinada em 
-20 dB/década, com um valor de constante negativo na curva de fase do diagrama. As Figuras 8.3 e 
8.4 ilustram os Diagramas de Bode para os dois casos. 


A potenciação do termo (jw)*! (por exemplo, de um número n, tal que o fator se torne (jw)*”, 
conduzirá a uma inclinação, respectivamente, de +20ndB/década nas curvas de módulo, que igual- 
mente se tornarão +20nlog(w). 


Módulo (dB) 
I 
N 
o 


0,001 0,01 0,1 1 10 


Fase (graus) 


0,1 


Log (frequéncia) 


Figura 8.3 - Diagramas de Bode para um fator jo. 
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Módulo (dB) 


Fase (graus) 


0 
-25 
— 50 
-75 
— 100 
— 125 
— 150 
- 175 


0,001 0,01 0,1 1 10 


Log (frequéncia) 


Figura 8.4 - Diagramas de Bode para um fator (jw»)”. 


A Figura 8.5 ilustra, de forma exemplificada, as curvas de módulo e fase para o termo (jo), 
em que se observa claramente uma inclinação de -60 dB/década para a curva de módulo. 


Projeto de Controladores no Domínio da Frequência - O Diagrama de Bode 


El 
3 
3 3 
2 = 
E 3 
= 
W 
`O 
0,001 0,01 0,1 1 10 
175 
150 HHHH 
7 125 
a 
5 100 
2 75 
lu) 
LL 
50 + HHHH 
25 
0 A+ 
0,001 0,01 0,1 1 10 
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Figura 8.5 - Diagramas de Bode para um fator (ja). 


+1 


8.2.3 Fatores de primeira ordem (1+ joT) 


O uso de fatores de primeira ordem (1+ joT)*! produz para as curvas de módulo do Diagrama 
1 
de Bode, respectivamente, os termos recíprocos 220198 7 =+20logvl+w?T? dB. 
+]0 


É interessante observar que tais termos possuem aproximações para frequências 
muito baixas ou muito altas. Quando a frequência se aproxima de zero, ou seja, quando 
w<1/T, o termo de primeira ordem em dB pode ser escrito para a curva de módulo como 
+20log V1+w?T? =+20log1=0 dB. Do mesmo modo, para frequências muito elevadas, ou seja, 
quando œ >1/T , o termo de primeira ordem em dB pode ser escrito para a curva de módulo como 
+20log V1+@?T? = +20logoT dB. Assim, a representação logarítmica da curva de resposta em fre- 
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quéncia do fator de primeira ordem (1+ jwT)*! pode ser aproximada por duas retas assíntotas, a 
primeira delas constante em 0 dB e a segunda com inclinação de + 20dB/década, respectivamente. 
A frequência em que o gráfico “abandona” as assíntotas é determinada “frequência de canto”, sendo 
iguala @ = T rad/s para o caso em questão. 


Para as curvas de fase referentes ao termo (1+ joT)*!, o valor da curva é representado pelos 
termos +tg !wT, respectivamente, que possuem igualmente valores peculiares nas frequências limí- 
trofes. De fato, para œ = 0 rad/s, +tg OT = 0; para a frequência de canto, +tg WT = 45º; por fim, 
para frequências muito elevadas, +tg !wT = 90º. Isso permite avaliar o comportamento das curvas 


de fase. As Figuras 8.6 e 8.7 ilustram o comportamento das curvas de módulo e fase para fatores de 
primeira ordem. 


T 
2 
2 
| 
3 
© 
> 
— 100 
T —120 
> 
fo) 
2 
g — 140 
lo] 
LL 
— 160 


Log (frequência) 


Figura 8.6 - Diagramas de Bode para um fator (1+ jwT), com T = 1s. 
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Módulo (dB) 


— 60 


Fase (graus) 
| 
8 


— 80 


aaa dra ab ddr dal 


0,1 1 10 
Log (frequéncia) 


Figura 8.7 - Diagramas de Bode para um fator (1+ jo T)”!, com T = 1s. 


8.2.4 Fatores quadráticos [1+22j0/0,)+G0/0,2]” 


Os termos quadráticos E +24(jo/@,)+(j@/ 0) | podem ser representados, em termos de dB, 
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2 
de módulo assume o valor de +20log > =+40 log dB, com inclinação de +40dB/década, res- 
0) 


n n 
pectivamente. É interessante observar que as assíntotas se cruzam em ©=0 , sendo esta, assim, a fre- 


quência de canto associada ao fator quadrático. As Figuras 8.8 e 8.9 ilustram o comportamento das 


curvas de módulo e fase para fatores quadráticos. 


60 t-t 


50 LIA 


40 


Módulo (dB) 


30 


20 


Fase (graus) 


ETICA E pit LJ! 


1 10 


Log (frequência) 


Figura 8.8 - Diagramas de Bode para um fator 1+2£(j0/0,)+(j0/0,)?,com q, = V6rad/s e t= 2, 


EN 
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Figura 8.9 - Diagramas de Bode para um fator [1+26Go/0,)+ (¡o/0, ? T , com 0, = V6rad/s e = 


sl- 


Fique de olho! 


Sistemas que não possuem polos nem zeros no semiplano lateral direito do plano complexo, ditos sistemas de fase 
mínima, têm função de transferência determinada de forma unívoca somente a partir da curva de módulo do Diagrama 
de Bode. Já sistemas que possuem polos e/ou zeros no semipleno lateral direito, ditos sistemas de fase não mínima, pos- 
suem função de transferência dependente da curva de fase do Diagrama de Bode. 
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~ Vamos recapitular? 


Foram apresentados, neste capítulo, conceitos referentes às técnicas de controle de processos basea- 
das no domínio da frequência, em especial, a aplicação do Diagrama de Bode para análise da resposta de 
sistemas. 


O capítulo abordou, entre outros elementos, a justificativa ao emprego do domínio da frequência 
como procedimento de controle alternativo às técnicas puramente temporais, enunciando um conjunto 
de passos para tal, de modo que se faça bom uso da funcionalidade da análise na frequência. 


O Agora é com você! 


1) Esboce os gráficos de módulo e fase para as funções de transferência em malha aber- 
ta apresentadas a seguir. Depois, determine também a resposta de cada sistema a 
entradas impulso unitário e degrau unitário, identificando-os como de fase mínima 
ou de fase não mínima. 


2 
a) GS) == 
3 
b) G E. 
(8) s+3 
1 
c) E 
d) Ge(s)=2 
e) Gals)j= l 
C 2 +4s+3 
4 
à Gya 
4 
g) => 
5 
i a T 


: 2) Para um sistema de segunda ordem em malha fechada, o que se pode afirmar quanto 
A à resposta em módulo, à medida que se varia o fator de amortecimento É ? Apresen- 
te um exemplo numérico e faça considerações a respeito disso. 
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88 


3) 


4) 


Obtenha o Diagrama de Bode da função de transferência da malha aberta 


K 
G(s)H(s) =——————.. A partir do diagrama traçado, verifique se o sistema em 
s(s+1)(3s+2) 


malha fechada é estável com K = 3. Justifique sua resposta. 

Considere um sistema com realimentação unitária cuja função de transferência de 
malha aberta seja dada por G(s)H(s) = ET Obtenha a resposta em regime per- 
manente de tal sistema quando o mesmo for submetido aos sinais de entrada: 

a) x(t)=sen(3t)u(t) 

b) x(t)=2sen(t)u(t) 
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Projeto de 
Controladores pelas 
Técnicas de 
Ziegler-Nichols 


Para começar 


Este capítulo tem por objetivo apresentar os conceitos de projeto de controladores com base na 
técnica de Ziegler-Nichols. Essa técnica corresponde a uma modificação da sintonia de controladores 
PID, por meio de uma heurística conhecida. 


As informações contidas neste capítulo trazem um complemento ao projeto de controladores PID, 
fornecendo um novo modo de conduzir os estudos apresentados no capítulo anterior, de forma que se 
aperfeiçoem as regras de projeto já enunciadas. 


9.1 A técnica de Ziegler-Nichols 


A técnica (ou método) de Ziegler-Nichols corresponde a um método heurístico de sintonia 
de um controlador PID. Tal procedimento foi desenvolvido por John G. Ziegler e Nathaniel B. 
Nichols, em 1942, sendo até hoje visto como clássico na Teoria de Controle, mas ainda utilizado 
em sua essência (ainda que com algumas ligeiras modificações, baseadas em ensaios práticos, que 
aperfeiçoaram as fórmulas propostas por Ziegler e Nichols). Os dois métodos básicos de ajuste 
de Ziegler-Nichols (método da resposta ao degrau unitário e método da realimentação por relé) 
buscam a obtenção de uma mesma resposta preestabelecida para o sistema em malha fechada, dis- 
tinguindo-se no que se refere à natureza da informação característica do processo abrangido por 
cada tipo de resposta. Uma vez obtidas as informações requeridas para cada método, faz-se uso de 
fórmulas simples para o cálculo de ganhos do controlador projetado. 


Amplie seus conhecimentos 


John G. Ziegler e Nathaniel B. Nichols estudaram, por volta da década de 1940, métodos de sintonia de controladores 
PID ainda antes do advento dos computadores pessoais. Embora não sejam, de fato, os formuladores do controle PID, 
são considerados, muitas vezes, os “pais” da teoria. 


A justificativa para o uso de técnicas empíricas como as de Ziegler-Nichols é evidente. Quando 
o projetista tem à sua disposição a descrição (modelo) matemática da planta que deseja controlar, 
várias técnicas, como as apresentadas nos capítulos anteriores, podem ser aplicadas, com o intuito de 
se determinar os parâmetros do controlador, para se obter as especificações em regimes transitório 
e permanente em malha fechada (tal como bem-definido na Figura 9.1). 


R(s) C(s) 
À E E 
PS a s ts +1 


Figura 9.1 - Atuação de um controlador PID sobre uma planta bem definida. 


Todavia, em alguns casos, a planta a ser controlada é de difícil caracterização, requerendo do 
projetista uma aproximação experimental para a sintonia do controlador projetado (conforme a 
Figura 9.2). 


C(s) 


Figura 9.2 - Atuação de um controlador PID sobre uma planta não definida. 


Na sintonia de um controlador PID, os métodos de Ziegler-Nichols buscam ajustar os pará- 
metros K,, T; e Ty a partir da observação da resposta a uma entrada conhecida (no método a ser 
abordado a seguir, um sinal degrau unitário) ou no ajuste do valor de ganho proporcional que asse- 
gure estabilidade marginal, para o caso de se utilizar somente um controlador proporcional isolado. 
O elenco de regras de Ziegler-Nichols, aplicável, sobretudo, na presença de sistemas a serem con- 
trolados que não apresentem uma descrição matemática precisa (muito embora possam também 
ser utilizados quando tais sistemas são bem-descritos), estabelece em si valores de K,, T; e tį que 
assegurem estabilidade à planta a ser controlada. Entretanto, parâmetros como sobressinal máximo 
podem atingir valores indesejáveis, de modo que os valores sugeridos de parâmetros nas técnicas 
de Ziegler-Nichols podem servir como pontos de partida para o ajuste fino de tais parâmetros, nem 
sempre adequados à primeira iteração. 


As seções seguintes descrevem cada método de Ziegler-Nichols, apresentando seus procedi- 
mentos e trazendo também exemplos de cada caso, para melhor compreensão dos passos de projeto. 
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9.2 Método da resposta ao degrau unitário 


O método da resposta ao degrau unitário (também referenciado como método do domínio 
temporal) exige o conhecimento de duas grandezas que caracterizam a resposta ao degrau de uma 
planta. Esse método consiste em se obter experimentalmente a resposta da planta a ser controlada, 
a uma entrada do tipo degrau unitário, como ilustrado na Figura 9.3. A não ser pela presença de 
polos na origem ou de polos complexos conjugados dominantes, a resposta ao degrau unitário terá a 
conformação de uma curva em S, sendo tal curva ajustada a partir de dois parâmetros: o atraso apa- 
rente L e a constante de tempo T. Tais parâmetros são obtidos pelo traçado de uma reta tangente à 
curva de resposta em seu ponto de inflexão, que cruza a reta c(t) = K, ao que os parâmetros L e T são 
determinados, conforme ilustrado na Figura 9.3. 


c(t) $ 
K 
> 
0 t 
— L T 


Figura 9.3 - Curva em S para o método de Ziegler-Nichols de resposta ao degrau unitário. 


A função de transferência “() pode ser definida como um sistema de primeira ordem dotado 
U(s) 
de atraso de transporte, conforme a Equação 9.1. 


C(s) Ke 
U(s) Ts+1 


(9.1) 


A Tabela 9.1 apresenta as regras de sintonia de Ziegler-Nichols, para o caso de análise da res- 
posta ao degrau da planta a ser controlada. 
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Tabela 9.1 - Regras de sintonia de Ziegler-Nichols, para o caso de análise da resposta ao degrau 


Segundo a Tabela 9.1, o controlador PID sintonizado de acordo com as regras apresentadas 
pode ser expresso, em sua planta, por: 


1 
Ge-(sS)=K,| 1+7,s+— 
c(s) o d 2) 


i 


T 1 
=1,2—| 1+0,5Ls + — 
A 21 (9.2) 


9.3 Método do período crítico 


O método do período crítico (também referenciado como método da sintonia por relé) exige o 
conhecimento de duas grandezas características da resposta em frequência da planta analisadas. Neste 
método, definem-se arbitrariamente 1, =% et, = 0, ajustando-se o valor do ganho proporcional K, 
(de modo que o controlador PID passe a se comportar como um controlador proporcional puro, con- 
forme ilustrado na Figura 9.4), até se obter o valor crítico Kop (para o qual começa a surgir oscilação, 
com período crítico T op), de maneira que a função de transferência do controlador PID se reduza a: 


1 
Go(s)=Kp, | 1L+t;s+— 
els o a 2) 


Jl 


1 
=0,6K ¿q | 1+0,125T pS + — (9.3) 
cal cr, T= y 
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R(s) C(s) 
+ Ko + planta > 


Figura 9.4 - Estratégia de controle para o método de Ziegler-Nichols do período crítico. 


Deve-se ressaltar que o valor de Kop pode ser obtido, caso o sistema a ser controlado possua 
planta matematicamente modelada, a partir da análise do Lugar das Raízes, bem como o período 
crítico de oscilação (Top). Tais valores são facilmente identificáveis nos gráficos de LGR, a partir das 
intersecções dos traçados com o eixo das frequências. 


A Tabela 9.2 apresenta as regras de sintonia de Ziegler-Nichols, para o caso de análise de 
ganho e período críticos. 


Tabela 9.2 - Regras de sintonia de Ziegler-Nichols, para o caso de análise de ganho e período críticos 


M RE E 


Tipo de controlador 


Exercício Resolvido 


1) Para o sistema de controle apresentado a seguir, aplique a regra de Ziegler-Nichols para 
projetar um controlador PID, de modo que o sobressinal máximo não exceda 25% do valor 
de regime permanente. 


R(s) A C(s) 
+ Gs) — — 
s(s + 2)(s + 3) 


Figura 9.5 - Exemplo de sistema de controle. 


Solução 


O controlador PID é expresso na forma: 


Il 
G (s)=K | l+t;s+— |. 
c(9) o d 2) 
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Para determinação dos parâmetros K,, Ôg e Ô, faz-se a escolha do segundo método de 
Ziegler-Nichols, já que a planta possui polos na origem. Fazendo ô} = 0 e ô = «o, a função 
em transferência de malha fechada se torna: 


d 


Oel Kp 
R(s) s(s+2)(s+3)+Kp 


C(s) _ Kp 
R(s) sê +582+65+K, 


Utilizando o Critério de Routh-Horwitz, é possível a determinação do valor de K, que gera 
oscilações não divergentes. 


6 
5o B 
30-K, 
5 
Kp 


O valor K, = 30 conduz a um elemento nulo na primeira coluna da tabela de Routh, de modo 
que esse é o valor crítico de ganho Kop: Isso conduz ao polinômio característico sî+5s%+6s+30. 
O valor da frequência das oscilações para o caso crítico (estabilidade marginal) é obtido pela 
substituição s = jœ no polinômio característico igualado a zero, o que resulta em: 

(jo) +50)? +6(j0)+30=0.. @= V6 =2,45. 


A partir desse valor de frequência, pode-se determinar o valor do período crítico de oscilação: 


Tu = 2 sta 


De posse dos valores de Kop e de Top obtém-se: 
Kp =0,6K ¿q = 18 
ta =0,125T¿g =0,321s. 


T =n = 1.2808 
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Isso conduz, passo a passo, a: 


1 
Go[(s)=K,|1+7,s+— 
si o É = 


1 
Go(s)=18| 1+0,321s+ 
c6) E) 


jl 2 +1,282s+1 
acos 603s2 +1,282s ) 


1,282s 

2 +0,8s +0,624 

A SEa : 
S 


O controlador PID possui 1 polo na origem e 2 zeros em s = -0,4 + j0,681. A função de trans- 
ferência em malha fechada é, portanto: 


1 
co Cras 
R(s) i 1 
ASE 


8 s2 +0,85+0,624 
C(s) É st + 5s3 + 68? 
R(s) 2 
(s) 148,755 s2 +0,85+0,624 
st 4583 + 68? 
Co 8,755(s2 + 0,85 + 0,624) 


R(s) s!+5s3+652+8,755(82 +0,85+0,624) 


C(s) _ 8,7558? + 7,004s + 5,463 
R(s) st+583 +14,755s2 +7,004s + 5,463 
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A resposta de tal função de transferência ao degrau unitário é ilustrada na Figura 9.6. 


20 30 40 50 60 70 80 90 100 


10 


t 


Figura 9.6 - Resposta ao degrau unitário para o sistema analisado. 
Uma vez que as especificações de sobressinal são atendidas (seu valor é aproximadamen- 


te 23,7% superior ao valor de regime), poder-se-ia tornar o sistema mais veloz alterando-se a 
posição dos zeros do controlador PID. Todavia, sendo o sobressinal a especificação 


única a ser 


z 


atendida, o controlador projetado náo necessita de ajustes finos em seus parámetros. 
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Fique de olho! 


O ajuste fino dos parâmetros após os valores inicias sugeridos pelo método de Ziegler-Nichols não é incomum. Todavia, 
há de se ter cuidado com os valores que se tomem para Kp, 1, € Ty, posto que não somente o sobressinal é afetado, mas 
também os valores de tempo de subida e de acomodação. Não é também proibitivo se alterar a localização de zeros no 
controlador PID projetado, atentando-se para o cuidado de não se descumprirem as especificações requeridas no projeto. 


( Vamos recapitular? 


Foram apresentados, neste capítulo, conceitos referentes aos métodos de Ziegler-Nichols para sin- 
tonia fina de controladores PID. Fazendo uso de regras empíricas de fácil manipulação, o emprego de 
tais técnicas permite que sistemas de difícil modelagem matemática possam ter especificações de desem- 
penho atendidas, na medida em que se faz uso de tabelas de aplicação imediata. 


O capítulo abordou, entre outros elementos, os dois métodos de aplicação de Ziegler-Nichols 
(método da resposta ao degrau unitário e método do período crítico), destacando em que situações 


são aplicáveis. 
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O Agora é com você! 


E 


Para os sistemas apresentados a seguir, projetar controladores PID G(s) baseados 
nas regras de Ziegler-Nichols, de modo que se atenda à especificação de sobressinal 
não superior a 20%. Se preciso for, faça ajustes finos nos parâmetros calculados. Uti- 
lize o método de Ziegler-Nichols (resposta ao degrau ou período crítico) que mais se 
aplique à configuração de polos da planta apresentada no problema de controle. 


R(s) 4 C(s) 
+ >| G¿(s) >» ———— m= 
s(s + 4)(s + 6) 
a) k 
R(s) 4 C(s) 
- Gas) he ———— HH 
b) (s + 3)(s + 7) 


R(s) 4 C(s) 
i Ki m? 
c) 
R 
À (s) A 
d) si + 252 + 3s +4 
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Análise de 
Controladores no 
Espaço de Estados 


Para começar 


Este capítulo tem por objetivo apresentar os preceitos de análise de sistemas de controle no espaço 
de estado, tópico bem-fundamentado na Teoria de Controle Moderno, sem adentrar de forma especí- 
fica o projeto de sistemas baseado em tal abordagem. Diferentemente de técnicas clássicas, que buscam 
relacionar saídas e entradas de sistemas, o controle moderno parte do pressuposto de que um sistema 
pode ser representado por uma única equação matricial de primeira ordem, agrupando todas as equa- 
ções diferenciais relacionadas ao sistema em questão e fornecendo um modo compacto de se modelar 
e analisar sistemas com múltiplas entradas e saídas. O projeto de controladores no espaço de estados é 
brevemente referenciado ao final do capítulo, na sugestão de técnica de alocação de polos, merecendo 
maior atenção em outras fontes bibliográficas que descrevem a teoria específica de Controle Moderno. 


10.1 Sistemas em espaço de estados 


Na representação de sistemas em espaço de estados, toma-se primeiramente como base uma 
equação que estabeleça relação entre saídas e entradas correspondentes a tais sistemas. Seja um sis- 
tema descrito pela equação diferencial 


@) (iD (n2) MM. ay (n2) 
yta y +8, y +..+a, ¡y+a,y=b,u+b, u +b, u +..+b, u+bu, (10.1) 


sendo y o sinal de saída e u o sinal de entrada. Tal equação pode ser representada no domínio da 
transformada de Laplace por: 


Ho dy +b A, eb, SD, (10.2) 
U(s) s+ a duda a 


A notação estabelecida na Equação 10.2 pode ser então representada por diferentes formas no 
espaço de estados, listadas a seguir. 


10.2 Representação de funções de transferência 
de sistemas no espaço de estados 


10.2.1 Forma canônica controlável 


A Equação 10.2 é representada na forma canônica controlável da seguinte forma: 


ži 0 0 o || x | fo 
Xy 0 0 1 o | x | jo 
p la] E A E (10.3) 
E 0 0 0 É lx, , | [0 
Xn “a, “aa “âno “MJ * 1 
Xx 
X5 
y =[b, —a,bo]D,-,-2, ¡bp | lbi aby ]| : |+bou (10.4) 
Xn-1 
x 


As representações utilizadas nas Equações 10.3 (equação de estados de entrada) e 10.4 (equa- 
ção de estados de saída) são úteis no projeto de sistemas de controle no espaço de estados baseado 
na alocação de polos, como brevemente enunciado neste capítulo. 


10.2.2 Forma canônica observável 


A Equação 10.2 é representada na forma canônica observável por meio das seguintes equações 
de estado: 
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X 000 =á Ži b, —3,B, 
Xy 0 Aral b,,¡ 72p-1Do 
o lse É: -a2 || © |+ : u (10.5) 
ža- 0 : ri b, -azb, 
Ža =a Xn b, ayb, 
xi 
X2 
y=[0 0 + 01] : |+bu (10.6) 
Xn- 
É 


Na Equação de entrada 10.5, a matriz de estados multiplicativa do vetor de entrada x é a 
matriz transposta da mesma utilizada na Equação de entrada 10.3. 


10.2.3 Forma canônica diagonal 


A forma canônica diagonal assume que os polos da função de transferência expressa na Equa- 
ção 10.2 são todos distintos, no que se reescreve a Equação 10.2 como: 


Y(s) _ bos" +bys" +... +b ,5+b, 
U(s) (s+p,Ms + p,)...(s+p,) 


(10.7) 


E A TR Cn 
U(s) S+P, S+P, S+Pn 


A Equação 10.7 é representada na forma canônica diagonal por meio das seguintes equações 
de estado: 


X -Pi X; 1 
X, —P, X, 1 
é |= É” 3 |+: fu (10.8) 
Ža- Xai 1 
X | =p. || E, 1 
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yla ES am G A : +bou (10.9) 


10.2.4 Forma canônica de Jordan 


A forma canônica diagonal, quando o sistema representado em questão apresenta raízes múlti- 
plas, necessita de uma pequena modificação, caso em que é denominada forma canônica de Jordan. 
Consideremos agora a Equação 10.7 como: 


Y(s) + sb toh E 
U(s) C+E CtP) Panna) Fpa) 
(10.10) 
Te r A E qa. 
U(s) (s+p,)" (asp e S+Pi S+Pm+ Ss+P, 


Amplie seus conhecimentos 


Camille Jordan 


O matemático Camille Jordan (1838-1922) iniciou seus 
estudos em matemática aplicada em 1869, motivado pelo 
interesse em estruturas cristalinas. Particularmente, Jordan 
tinha gosto peculiar por aprofundar seu conhecimento em 
transformações euclidianas tridimensionais. Intensificando 
suas pesquisas em Teoria de Grupos, Jordan tornou-se 
um nome de destaque na formulação de subgrupos finitos 


Wikimedia/Commons 


de um grupo linear de matrizes quadradas no espaço dos 
números complexos. Seu nome é lembrado em definitivo, 
porém, quando se associa a ele a prova de que uma curva 
fechada divide um plano em exatas duas regiões (o cha- 
mado teorema da curva de Jordan). 


Em Controle, Jordan trouxe sua contribuição em 1870, 

quando publicou um teorema sobre matrizes quadradas que 

provê a base de comprovação de outros teoremas simples e Figura 10.1 - O matemático Camille Jordan, 
que permite a representação de vários sistemas lineares na que grande contribuição trouxe à análise 
forma de equações matriciais, independentemente da multi- de sistemas no espaço de estados. 
plicidade de polos. 


As equações de entrada e de saída para a forma canônica diagonal, representadas, respectiva- 
mente, pelas Equações 10.8 e 10.9, tornam-se: 
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X; | -p 1 
X 0 -pP 
X; 0 
Xa 

21 10 


Y =| C2 


1 
“Pij 0 
0 -Pma 
0 0 
Cn 6, ] 


u (10.11) 


+. oso 


(10.12) 


Observa-se que as Equações 10.10 e 10.12 permanecem idênticas, posto que a forma canónica 
de Jordan provoca, substancialmente, uma alteração na representação dos polos do sistema. A téc- 
nica nada mais é do que uma forma alternativa de se representar uma matriz ou operador linear por 
meio de outra matriz semelhante à original, sendo quase uma matriz diagonal. 


A contribuição efetiva da forma canônica de Jordan reside no fato de que seu uso permite que 
se extraiam informações sobre a transformação aplicada a sinais de entrada com facilidade, lem- 
brando que é necessário que se conheça a multiplicidade algébrica de cada um dos autovalores asso- 
ciados para que seja efetiva a utilização da forma canônica de Jordan. 


1) Para o sistema dado por Y(S) _ 


U(s) s2+7s+12 


s+2 


nas formas canônica controlável, observável e diagonal. 


Utilizando as equações anteriores, podemos escrever: 


a) Forma canônica controlável 


A equação de estado de entrada dessa forma canônica é dada por: 


É o 1 0 
E, o 0 1 

É 0 0 0 
a | an an- za n-2 
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, obtenha sua representação no espaço de estados, 


Observando que a, = 7 e a, = 12, tem-se: 


H 


A equação de estados de saída é definida por: 


0 


š, | |-12 


y=[b, -a,bo |b -a, Do |---|b; FE 


y= [b, -abọ b; ado] 


abo | 


2 
+ bu 
x 
2 


Considerando que se tem os valores obtidos dos coeficientes, a, =7, a, = 12,b,=0,b, =1,b,=2, 


escreve-se: 


zi 


y =| -24 ais] 


x7 


b) Forma canônica observável 


A equação de estado de entrada dessa forma canônica é 
X; 0 E X; 
x, 1 A S 
5 = ES : 
Xa- Xa 
Xn aj n 
o que conduz a 
X =|| ss 
Xx, OA ES 


| 


dada por: 
| ba —anDo 
Da- 7 an-1P0 
+ E u, 
b, —a,bo 
b; =ajb, 
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A equação de estados de saída é obtida por: 


y=[0 0: 01] : |+bu 


c) Forma canônica diagonal 


A equação de estado de entrada dessa forma canônica é dada por: 


ES =P 2 
x7 =P X 
2 = +|: 
Sl Snl 
Xa "Pr Xa 


E S 


Igualmente, a equação de estado de saída é obtida de: 


yie E do e Cal 2 Ebu 
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10.3 O método da alocação de polos 


O método da alocação de polos é, de certo modo, similar ao método do LGR, em que se defi- 
nem as posições exatas de todos os polos de malha fechada nas posições desejadas. No método de 
alocação de polos, admite-se que, se o sistema considerado for de estado completamente controlável 
(pode partir de um estado inicial qualquer a outro final dentro de um intervalo de tempo desejado), 
será possível alocar os polos de malha fechada em qualquer posição, com o uso de realimentação 
apropriada. 


Nessa técnica, determinam-se os polos de malha fechada requeridos, tomando-se como base 
as especificações dadas de resposta temporal e/ou em frequência, por exemplo, tempos de subida e 
coeficiente de amortecimento, bem como especificações de regime permanente, como o tempo de 
acomodação. Usualmente, o projeto de um sistema de controle toma como base o conhecimento 
prévio da localização dos polos de malha fechada, visando atingir-se um dado coeficiente de amor- 
tecimento e uma frequência natural preestabelecida - na condição de que polos não dominantes não 
exerçam algum efeito considerável sobre a saída do sistema. A técnica de alocação de polos admite 
que todos os polos de malha fechada, e não apenas os dominantes, sejam previamente alocados, 
desde que a controlabilidade total se faça presente. 


O projeto específico de controladores baseados no espaço de estados extrapola o escopo desta 
obra, mas um estudo apurado do tema pode ser facilmente conduzido, na condição de que o leitor se 
aprofunde no estudo da álgebra linear. 


Fique de olho! 


O projeto de controladores no espaço de estados é particularmente útil em problemas nos quais um bom modelo mate- 
mático do sistema a ser controlado esteja disponível e em que a dimensão dos estados não seja muito elevada. Nessas 
condições, a técnica se apresenta como de mais rápida implementação do que as técnicas apresentadas na teoria clás- 
sica de controle, exigindo, entretanto, o conhecimento de mais parâmetros do que os usuais. 
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——( Vamos recapitular? 


Foram apresentados, neste capítulo, conceitos referentes à análise de sistemas de controle no 
espaço de estados, sem adentrar especificamente o modelamento de controladores. 


O capítulo abordou, entre outros elementos, as formas canônicas controlável, observável e diago- 
nal, exemplificando a transformação de equações diferenciais em equações matriciais. 


O: Agora é com você! 


Para os sistemas representados a seguir, obtenha sua representação no espaço de 
estados, nas formas canónica controlável, observável e diagonal (ou de Jordan, se 
possível). Determine também os valores de coeficiente de amortecimento de tais 
sistemas, escolhendo um deles para o qual SS > 50%, e projete um controlador PID 
conforme a técnica apresentada no capítulo anterior, para que, em malha fechada, o 
novo valor de SS náo exceda 30%. 


Y (s) s+4 
O) == 
U(s) s2+3s+4 
b Y(s) _ 8? +25+2 
) U(s) s?+3s+2 
Y(s) _  s?+4s+2 
e) U(s) s+482+75+12 
Y(s s 
d) Y(s) _ 
U(s) s2+12 
Y(s) s2+45+2 
O) === 


U(s) s3+6s+1 
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Y(s) _ s2+4 


0 U(s) s3+382+3s+1 
Yo 2 

O To) sA 

h) Y(s)_  s+2 


U(s) s?-4s+2 


Sob que condições a forma canônica de Jordan se apresenta como mais vantajosa em 
relação à forma diagonal simples? Exemplifique com situações práticas tal escolha. 
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Aplicações de 
Sistemas de 
Controle em 

Processos 
Industriais 


Para começar 


Este capítulo tem por objetivo apresentar nichos de aplicação de sistemas de controle, focando 
majoritariamente em processos industriais cujas estratégias de controle são de usual aplicação. 


Os tópicos aqui apresentados não fazem ênfase a uma ou outra técnica de controle descrita nos 
capítulos anteriores, enfatizando a gama de possibilidades de uso de sistemas de controle a diversas plan- 
tas e seus elementos. Dá-se ênfase a dois elementos já descritos no capítulo inicial desta obra: atuadores 
e sensores, descrevendo-se como estratégias de controle que se aplicam à operação de tais dispositivos 
em uma planta. Não se busca aqui trazer uma discussão aprofundada sobre a natureza de tais elementos 
- posto ser essa uma tarefa mais afeita aos preceitos da Robótica -, mas introduzir conexões entre seções 
físicas de plantas e sistemas de controle a ela associados. 


11.1 Introdução às aplicações de Sistemas de Controle 


As técnicas de projeto de Sistemas de Controle descritas nos capítulos anteriores se prestam a 
uma variedade ampla de aplicações, notadamente em processos que envolvem um elevado número 
de etapas, caso em que plantas industriais se encaixam. De um modo objetivo, este capítulo traz uma 
correlação entre as técnicas de controle apresentadas e as mais usuais aplicações industriais, em ter- 
mos dos elementos mecânicos atuadores (componentes que realizam a conversão de energia elétrica, 
hidráulica ou pneumática em energia mecânica) e sensores (transdutores de grandezas físicas em 
sinais elétricos correspondentes). 


11.2 Elementos de atuação (atuadores) 


Atuadores são elementos encontrados em plantas representativas de processos caracteriza- 
dos substancialmente por produzir algum movimento, sempre comandados por sinais de natureza 
manual, mecânica (nas modalidades hidráulica ou pneumática) ou elétrica. Na indústria, um atua- 
dor é um mecanismo utilizado de forma comum na fabricação de máquinas e equipamentos que se 
prestam a iniciar ou finalizar funções comandadas por válvulas. 


De modo a ressaltar o papel de sistemas de controle sobre o funcionamento de tais dispositi- 
vos, as próximas subseções trazem de forma objetiva os tipos de atuadores encontrados na indústria. 


11.2.1 Atuadores manuais 


Atuadores manuais correspondem a dispositivos que fazem uso de alavancas, engrenagens ou 
polias para facilitar o movimento, desacoplados de uma fonte externa de energia que produza tal 
movimento. É comum o uso de atuadores manuais em tarefas que exijam pouco esforço mecânico, 
geralmente atrelados a válvulas e outros equipamentos de pequenas dimensões. Todavia, não são 
recomendados atuadores manuais em ambientes hostis ao ser humano, como em ambientes tóxicos 
ou insalubres. A Figura 11.1 ilustra um atuador manual. 


Valveactuator/Wikimedia Commons 


Figura 11.1 - Atuador (cilindro) manual. 


A ação básica de controle de atuadores manuais é notadamente representada por um sistema 
em malha aberta, posto que a observação da saída não produz diretamente um erro a ser realimen- 
tado à entrada, de modo a controlar o sinal reinjetado no sistema. 
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11.2.2 Atuadores hidráulicos 


Atuadores hidráulicos são aqueles que fazem uso de fluidos para pressurizar pistões, tor- 
nando fácil a realização de dado esforço mecânico. Pelo fato de fluidos hidráulicos não poderem 
ser comprimidos, variáveis de processo como velocidade e potência são de mais difícil obtenção, 
permitindo-se, todavia, a sintonia fina de precisão no movimento controlado. Usualmente, atuado- 
res hidráulicos são representados pelo cilindro hidráulico, como mostra a Figura 11.2, cuja função 
primeira é a de transformar energia hidráulica em energia mecânica. O cilindro hidráulico é parte 
constituinte de diversos sistemas de manipulação mecânica e articulação, desempenhando a aplica- 
ção de uma força em um percurso linear, gerando o movimento de um pistão móvel. Tal força é 
constante ao longo do movimento do pistão (assim como ao produto da área da seção transversal do 
cilindro pela pressão da bomba). 


Cilindros hidráulicos possuem aplicação conjunta com elementos de maior dimensão, por 
exemplo, no movimento de pás de retroescavadeiras. São também empregados em outros setores em 
que se exijam forças elevadas, em detrimento de velocidade. 


A ação básica de controle de atuadores pneumáticos reside no fato de que se deva privilegiar 
precisão (ou seja, mínimo erro em regime estacionário), o que pressupõe uma ação de controle pro- 
porcional integral, não sendo a velocidade (ou seja, o tempo de subida é o mesmo de acomodação) 
privilegiada na análise. 


Rstom03/Wikimedia Commons 


Figura 11.2 - Atuador (cilindro) hidráulico. 
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11.2.3 Atuadores pneumáticos 


Atuadores pneumáticos são aqueles que utilizam ar para pressurizar os pistões, em vez de 
outro fluido. Como o ar pode ser comprimido, dispensam regulação prévia; além disso, são mais 
rápidos nas situações de partida e parada, o que os recomenda para operar em equipamentos mecá- 
nicos que requerem mudanças repentinas e bruscas de pressão para executar uma dada tarefa. A 
Figura 11.3 apresenta um cilindro pneumático que representa tal categoria de atuador. 


Mixabest/Wikimedia Commons 


Figura 11.3 - Atuador (cilindro) pneumático. 


11.2.4 Atuadores elétricos 


Atuadores elétricos são aqueles movidos por um motor, o qual fornece torque para ope- 
rar válvulas em um equipamento mecânico. São comumente utilizados quando o equipamento 
a realizar movimento requer válvulas multiturno. Uma vez que essas válvulas são utilizadas na 
presença de maquinário com frequente mudança de marcha em funcionamento, atuadores elé- 
tricos sofrem razoável demanda, além de requerer alimentação de backup para funcionamento 
ininterrupto. 


Eles são costumeiramente referidos como motores elétricos - embora a denominação não seja 
fisicamente inconsistente, prefere-se, aqui, distinguir tais elementos, conquanto motores provoquem 
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o movimento em si nos atuadores, em vez de representá-los inteiramente. Nesse sentido, motores 
podem ser subdivididos em: 


» 


» 


» 


» 


Motores de corrente alternada (AC): convertem energia elétrica em energia mecânica. 
Motores de corrente contínua (DC): fazem uso de uma fonte de tensão contínua. 


Motores de passo: comportam-se como motores de corrente contínua desprovidos de 
escovas e comutadores, dotados de altíssima precisão linear ou angular. 


Servomotores: consistem em motores de passo com pouca capacidade de realizar elevados 
esforços; são dotados de um sistema de controle interno que verifica a posição de entrada 
em relação à posição de saída. 


As Figuras 11.4 e 11.5 apresentam alguns tipos de atuadores elétricos. 
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Figura 11.4 - Atuador elétrico - motor de passo. 


José Luis Gálvez (Digigalos)/Wikimedia Commons 


Figura 11.5 - Atuador elétrico - servomotor. 


Fique de olho! 


Um servomotor é bastante empregado na otimização da operação de máquinas e equipamentos. 


Atualmente, dispõe de elementos de realimentação (malha fechada), em que cada servomotor é acompanhado por um 
servoconversor, o que permite, por exemplo, o controle do torque de um eixo de máquina para amplas faixas de velocida- 
de de rotação. 


Assim, servomotores encontram aplicação tanto em tarefas de pequeno porte quanto em tarefas mais pesadas na 


indústria, sobretudo quando for desejado o controle em malha fechada de variáveis de processo, como torque, posição 
ou velocidade. 


11.3 Elementos de sensoriamento (sensores) 


Sensores são elementos utilizados para medir e monitorar variáveis de processo, como tempe- 
ratura, vazão, nível, pressão, entre outras, tendo surgido quando da necessidade de automatização 
dos processos, visando a uma maior eficiência e eficácia dos processos. Podem ser fisicamente des- 
critos como transdutores que alteram sua constituição interna diante da ocorrência de um fenômeno 
físico externo, convertendo tais estímulos externos em sinais elétricos que atuam sobre a variável de 
processo causadora desses estímulos. Atualmente, sensores possuem uma extensa gama de aplica- 
ções, existindo de variadas formas e modelos. As subseções seguintes descrevem, de forma breve, 
quatro tipos de sensores, bastante empregados na indústria. 


11.3.1 Sensores de nível 


Sensores de nível são aqueles que se prestam a mensurar a altura do conteúdo de um reserva- 
tório qualquer, que armazene matéria líquida ou sólida. O sensoriamento de nível pode se dar por 
medição direta, cuja referência de medição é o limite superior da substância medida (por exemplo, 
pelo uso de réguas, gabaritos, visores ou boias), ou por medição indireta, cujo nível é determinado 
com base na medição de outras variáveis de processo, como pressão, radiação, empuxo etc. Veja na 
Figura 11.6 um exemplo de sensor de nível. 


Gillsensors/Wikimedia Commons 


Figura 11.6 - Sensor de nível. 
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11.3.2 Sensores de temperatura 


Sensores de temperatura são aqueles que se prestam a mensurar a temperatura de um dado 
meio ou matéria, convertendo-a em sinal elétrico. Usualmente, dividem-se em: 


» Sensores resistivos: resistências dependentes da temperatura. 
» Sensores termoelétricos ou termopares: valem-se do efeito Seebeck. 


» Sensores de infravermelho: captam radiação eletromagnética no comprimento de onda da 
radiação infravermelha, sendo tal faixa de frequência característica de emissão de calor. 


A Figura 11.7 ilustra um sensor de temperatura. 


Harke/Wikimedia Commons 


dae > 0 


Figura 11.7 - Sensor de temperatura - termopar tipo K. 


Amplie seus conhecimentos 


Thomas Johann Seebeck 


O Físico Thomas Johann Seebeck (1770-1831) foi o descobridor, 
em 1821, do efeito que leva seu nome, o qual descreve a produ- 
ção de uma diferença de potencial entre duas junções de condu- 
tores (ou semicondutores) de materiais diferentes, quando estas 
estão sujeitas a diferentes temperaturas. Juntamente com outro 
efeito físico — Peltier —, o efeito Seebeck descreve com exatidão o 
funcionamento de termopares. 
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Figura 11.8 - O físico Thomas Johann Seebeck. 


11.3.3 Sensores de pressáo 


Sensores de pressáo sáo aqueles que se prestam a converter uma pressáo medida em um 
sinal mecánico ou elétrico, constituídos por um elemento primário (que sofre ação - deformação 
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ou deflexão -, direta ou indireta, da pressão medida) e por um elemento secundário (que traduz 
variações de pressão em valores adequados a medição e controle). 


Podem ser classificados em: 

»  Manométricos: tomam como referência a pressão atmosférica. 

» Absolutos: tomam como referência a pressão do vácuo. 

» — Diferenciais: tomam como referência outro valor de pressão previamente conhecido. 


Em termos construtivos, sensores de pressão são usualmente encontrados como medidores de 
coluna líquida, medidores por deformação ou medidores eletroeletrônicos. 


A Figura 11.9 mostra um sensor de pressão manométrico. 


Zaphod/Wikimedia Commons 


Nro Oz 


Figura 11.9 - Sensor de pressão - manômetro. 


11.3.4 Sensores de vazão 


Sensores de vazão são aqueles que se prestam a medir a quantidade de material líquido ou 
gasoso que atravessa determinado ponto de uma tubulação durante certo intervalo de tempo. Os 
principais tipos de sensores de vazão se baseiam na medição indireta de pressão diferencial (por 
exemplo, placas de orifício, tubos de Venturi e tubos de Pitot), enquanto outros fazem uso de turbi- 
nas ou de campos magnéticos. A Figura 11.10 apresenta um sensor de vazão. 
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Figura 11.10 - Sensor de vazão - tubo de Pitot inserido na extremidade de uma aeronave. 


Fique de olho! 


Exemplo interessante do uso de sistemas de controle em plantas industriais é o da válvula redutora de pressão, instala- 
da em sistemas de tubulação conectados a tanques de armazenamento ou reservatórios de líquido ou gás. Essa válvula 
gerencia a pressão no sistema controlado por meio de um conjunto mola-diafragma, que abranda o fluxo de líquido no 
sistema, mantendo um nível preajustado de pressão. 


Válvulas redutoras de pressão de água são, muitas vezes, instaladas em sistemas residenciais de água, evitando o 
desperdício que se observaria caso toda a alta pressão de bombeamento de água das ruas fosse utilizada também 
nas casas. 


Vamos recapitular? 


Foram descritos neste capítulo elementos que servem de suporte para a aplicação de sistemas de 
controle em processos ou plantas industriais: os atuadores e os sensores. A operação de tais elementos 
segue as premissas de controle descritas nos capítulos anteriores, e as características de operação de 
cada um dos atuadores e sensores apresentados devem se adequar às especificações de desempenho 
de cada problema. 
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. 


Agora é com você! 


1) 


2) 


3) 


Aponte as diferenças básicas entre atuadores e sensores, descrevendo sua importân- 
cia em uma malha de controle completa. 


Explique o funcionamento de tubos de Pitot em aeronaves, destacando sua finalidade 
e os efeitos adversos de seu mau funcionamento. 


Você acha que o corpo humano é dotado de sensores e atuadores? Forneça exemplos 
para sua resposta. 
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Relações 


Trigonométricas 
Importantes 
1) tga = "2 10) ary ps 
cosa 
2) sen?a+cos?a=1 1) sen?a — 1700824 
3) seca= 
cosa 12) sen2a = 2senacosa 
4) coseca= Y 13) cos2a = cos? a -sen?a =1-—2sen?a = 2cos? a—1 
sena 
1 14) teo 2 
5) cotga=— I-tg?a 
tga 
6) sec?a=l+tg?a 15) font) = I-cosa 
2 2 
7) cosec? a =1+cotg%a 
a l+cosa 
16) |cos—|=,|—— 
8) sen(atb)=sena.cosbtsenb.cosb 2 2 
tgatt = 
9) te(atb)= 825180 gb 17) "e. cosa 
I-tgatgb 2 sena 


18) senacosb= 5 [senta +b) +sen(a —b)] 
19) senasenb = [costa —b)-cos(a +b)] 
20) cosacosb= S[eosta +b) + cos(a —b)] 
21) cosasenb = [senta +b) —sen(a —b)] 
22) sena-senb= Zsen( 257 )eos( 22) 


23) lsena=Itcos[ Z-a) 
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Arcos Elementares (Primeiro Quadrante) 


Funções 
Trigonométricas 
e Propriedades 
de Números 
Complexos 


Propriedades de Números Complexos 


7) 


8) 


9) 


A 
Im 


a+jb 


Re 


z=a+jb=rei? =r(cos0+ jsen0), r=va2 +b2, RA 2 =-1 
a 


r=|z]= Va? +b? 
z=a+jb, z*=a-jb 
z+z*=2. Refz) 

z.z* =|z? 


2125) =|21)-|z5| 


Vz =Vz[ cos 22n 
n 
z+z* 
Refzl=—— 
tz 3 
Z—2z* 
Im{z} = 
tz) 2j 


+ jse 


n 


0+2kr 
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Tabelas de Pares 
de Transformada 
de Laplace 


F(s) 


n! 
Ng-at = === 
E i d a 
(0) 
sen(wt 
CI | so 
S 
10 cos(wt 
OIEA 
(0) 
S 
1 
13 
s(s+a) 
1 1 
14 —— (erat etbt Ca a 
E Es ERR) (s+a)(s+b) 
1 S 
15 b -bt _ —at Et 
o EA á ae (s+a)(s+b) 


e pa Es 
7 [t Bi J ES GE 
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Propriedades 
da Transformada 
de Laplace 


L[Af(t)]= AF(s) 
L[ Af, (t) + Bf, (t) ] = AF (s) + BE,(s) 
qa. sE(s)—£(0,) 


EQ] o x 
c K F(s)-sf(0,)-f'(0,) 


c| fr (at |= = (condições iniciais nulas) 


d fr oa = 9 (condições iniciais nulas) 
n 


€ 


L| f(t) | =F(s+a) 


£[f(t-a)u(t-a)]= e-*F(s) 


c[errço]=(Dº Zro (n =1,2,3...) 


d2] =aF(as) 


a 


t 
MORAO -d j fof Ct- as =F (s)-E,(s) 
0 


lim f(t) = lim sF(s) 
to s>0 


